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1. Introduction. This paper provides the mathematical foundation for 
stochasticahy continuous affine processes on the cone of positive semidefi- 
nite symmetric d x d-matrices Sj" . These matrix- valued affine processes have 
arisen from a large and growing range of useful applications in finance, in- 
cluding multi-asset option pricing with stochastic volatility and correlation 
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structures, and fixed- income models with stochastically correlated risk fac- 
tors and default intensities. 

For illustration, let us consider a multi-variate stochastic volatility model 
consisting of a d-dimensional logarithmic price process with risk-neutral dy- 
namics 

(1.1) dYt = irl-lxf''^)dt + y^tdBt, Yo = y, 

and stochastic covariation process X = {Y, Y) , which is a proxy for the in- 
stantaneous covariance of the price returns. Here B denotes a standard d- 
dimensional Brownian motion, r the constant interest rate, 1 the vector 
whose entries are all equal to one and X'^^^^ the vector containing the diag- 
onal entries of X. 

The necessity to specify X as a process in S'j" such that it qualifies as 
covariation process is one of the mathematically interesting and demanding 
aspects of such models. Beyond that, the modeling of X must allow for 
enough flexibility in order to reflect the stylized facts of financial data and 
to adequately capture the dependence structure of the different assets. If 
these requirements are met, the model can be used as a basis for financial 
decision-making in the area of portfolio optimization, pricing of multi-asset 
options and hedging of correlation risk. 

The tractability of such a model crucially depends on the dynamics of X. 
A large part of the literature in the area of multivariate stochastic volatility 
modeling has proposed the following affine dynamics for X: 

dXt = {b + HXt + XtH^) dt + y^tdWt S + dW^ y^t + dJt, 

(1-2) 

Xo = xeSj, 

where b is some suitably chosen matrix in H, S some invertible matrices, 
W a standard d x d-matrix of Brownian motions possibly correlated with 
B, and J a pure jump process whose compensator is an affine function of 

The main reason for the analytic tractability of this model is that, under 
some technical conditions, the following affine transform formula holds: 

for appropriate arguments z G 5^ x iS^ and v G C^. The functions and ^ 
solve a system of nonlinear ordinary differential equations (ODEs), which 



■^This affine multi-variate stochastic volatility model generalizes the well-known one- 
dimensional models of Heston [27], for the diffusion case, or the Barndorff-Nielsen Shepard 
model [2], for the pure jump case. 
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are determined by the model parameters. Setting v = 0, 4>it,z) = —^{t,z,0) 
and xlj{t, z) = — ^(t, z, 0) and taking z = n G Sj", we arrive at 

(1.3) E.fe-T--^"^*)] = e-*(*'-)-Tr('/'{*.-)-), u G S+ . 

In this paper, we characterize the class of all stochastically continuous 
time-homogeneous Markov processes with the key property (1.3) — henceforth 
called affine processes — on Sj". Our main result shows that an affine pro- 
cess is necessarily a Feller process whose generator has affine coefficients 
in the state variables. The parameters of the generator satisfy some well- 
determined admissibility conditions, and are in a one-to-one relation with 
those of the corresponding ODEs for (f) and ^. Conversely, and more impor- 
tantly for applications, we show that for any admissible parameter set there 
exists a unique well-behaved affine process on . Furthermore, we prove 
that any stochastically continuous infinitely decomposable Markov process 
on S'^ is affine with zero diffusion, and vice versa. 

On the one hand, our findings extend the model class (1.2), since a 
more general drift and jumps are possible. Indeed, we allow for full gen- 
erality in 6, as long as b — {d — l)YP"ll G S^, for a general linear drift part 
-^(^) = Ylij^ijP^^ ^'^d for an inclusion of (infinite activity) jumps. This of 
course enables more flexibility in financial modeling. For example, due to 
the general linear drift part, the volatility of one asset can generally de- 
pend on the other ones, which is not possible for B{x) = Hx + xH~^ . On 
the other hand, we now know the exact assumptions under which affine pro- 
cesses on actually exist. Our characterization of affine processes on Sj" 
is thus exhaustive. Beyond that, the equivalence of infinitely decomposable 
Markov processes with state space Sj" and affine processes without diffusion 
is interesting in its own right. 

This paper complements Duffie, Filipovic and Schachermayer [16], who 
analyzed time-homogeneous affine processes on the state space W]^ x R".^ 
Matrix-valued affine processes seem to have been studied systematically for 
the first time in the literature by Bru [5, 6], who introduced the so called 
Wishart processes. These are generalizations of squares of matrix Ornstein- 
Uhlenbeck processes, that is, of the form (1.2) for J = and b = kYl^T,, for 
some real parameter k > d — 1. Note that A; > d — 1 is a stronger assumption 
than what we require on b and Bru [6] then establishes existence 

and uniqueness of a local^ S'J'-valued solution to (1.2) under the additional 
assumptions that Xq has distinct eigenvalues, —H G Sj", and that H and S 
commute (see [6], Theorem 2"). In the more special case where H = and 



''For the diffusion case see also [20] or [19], Chapter 10. Time-inhomogeneous affine 
processes on x R" have been explored in [18]. 
^Up to the first collision time of the eigenvalues. 
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k> d—1, Bru [6] shows global existence and uniqueness for (1.2) for any Xq 
with distinct eigenvalues (see [6], Theorem 2 and last part of Section 3).^ 
Bru's results concerning strong solutions have recently been extended to the 
case of matrix valued jump-diffusions; see [40]. 

Wishart processes have subsequently been introduced in the financial lit- 
erature by Gourieroux and Sufana [24, 25] and Gourieroux et al. [23]. Fi- 
nancial applications thereof have then been taken up and carried further by 
various authors, including Da Fonseca et al. [9-12] and Buraschi, Cieslak 
and Trojani [7, 8]. Grasselli and Tebaldi [26] give some general results on 
the solvability of the corresponding Riccati ODEs. Barndorff-Nielsen and 
Stelzer [3] provide a theory for a certain class of matrix- valued Levy driven 
Ornstein-Uhlenbeck processes of finite variation. Leippold and Trojani [38] 
introduce S'j'-valued affine jump diffusions and provide financial examples, 
including multi-variate option pricing, fixed-income models and dynamic 
portfolio choice. All of these models are contained in our framework. 

We want to point out that the full characterization of positive semidefi- 
nite matrix- valued affine processes needs a multitude of methods. In order 
to prove the fundamental property of regularity of affine processes another 
adaption of the famous analysis of Montgommery and Zippin is necessary, 
which has been worked out in [34] and [35] for the state space M™ x R". For 
the necessary conditions on drift, diffusion and jump parameters we need 
the theory of infinitely divisible distributions on Sj" . Most interestingly, the 
constant drift part b must satisfy a condition depending on the magnitude 
of the diffusion component (see Proposition 4.18), which is in accordance 
with the choice of the drift in Bru's work [6] on Wishart processes, as ex- 
plained above. This enigmatic additional condition on the drift b is derived 
by studying the process with respect to well chosen test functions, including 
in our case the determinant of the process. It is worth noting, as already 
visible in dimension one, that a naive application of classical geometric in- 
variance conditions does not bring the correct necessary result on the drift 
but a stronger one. Indeed, take a one-dimensional affine diffusion process 
X solving 

dXt = bdt+ ^/YtdWt. 

Then a back-of-the-envelope calculation would yield the Stratonovich drift at 
the boundary point x = of value b — j, leading to the necessary parameter 
restriction b> j, which is indeed too strong. It is well known that the correct 
parameter restriction is b>0. We see two reasons why geometric conditions 



^Actually, Bru [6] establishes existence and uniqueness of solutions also for k — 
l,...,d— 1. But these are degenerate solutions, as they are only defined on lower- 
dimensional subsets of the boundary of S'^ (see [6], Corollary 1). 
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on the drift cannot be applied: first, precisely at the boundary of our state 
spaces the diffusion coefficients are not Lipschitz continuous anymore, and, 
second, the boundary of the cone of positive semi-definite matrices is not a 
smooth submanifold but a more complicated object. 

For the sufficient direction refined methods from stochastic invariance the- 
ory are applied. Having established viability of a particular class of jump- 
diffusions on Sj", existence of affine processes on S'j" — under the necessary 
parameter conditions — is shown through the solution of a martingale prob- 
lem. Uniqueness follows by semigroup methods which need the theory of 
multi-dimensional Riccati equations. 

Summing up, we face two major problems in the analysis of positive ma- 
trix valued affine processes. First, the candidate stochastic differential equa- 
tions necessarily lead to volatility terms which are not Lipschitz continuous 
at the boundary of the state space. This makes every existence, uniqueness 
and invariance question delicate. Second, the jump behavior transversal to 
the boundary is of finite total variation. 

1.1. Program of the article. For affine processes on S'^ , results and proofs 
deviate in essential points from the theory on state spaces of the form x 
given in [16, 34], which is a consequence of the more involved geometry 
of this nonpolyhedral cone. The program of the paper as outlined below 
therefore includes a comparison with the approach in [16]. 

Section 2 contains the main definition and a summary of the results of 
this article. In Section 3, we then derive two main properties, namely the 
regularity of the process and the Feller property of the associated semigroup. 
The Feller property, in turn, is a simple consequence of an important posi- 
tivity result of the characteristic exponents which is proved in Lemma 
3.3. This lemma is further employed as a tool for the treatment of the gener- 
alized Riccati differential equations in Section 5.1 (see proof of Proposition 
5.3). The global existence and uniqueness of these equations is then used to 
show uniqueness of the martingale problem for affine processes (see proof of 
Proposition 5.9). 

In Section 4, we define a set of admissible parameters specifying the in- 
finitesimal generator of affine semigroups and prove the necessity of the 
parameter restrictions (see Proposition 4.9). 

The sufficient direction is then treated in Section 5. It is known that, for 
d>2, there exist continuous affine processes on which are — in contrast 
to those on the state space M!p x — not infinitely divisible (see Example 
2.8). The analysis of this paper reveals the failure of infinite divisibility as 
a consequence of the drift condition (see proof of Theorem 2.9). This has 
substantial influence on the approach chosen here to prove existence of affine 
processes associated with a given parameter set: Being in general hindered 
to recognize the solutions of the generalized Riccati differential equations as 
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cumulant generating functions of sub-stochastic measures, as done in [16], 
Section 7, we solve the martingale problem for the associated Levy type 
generator on S'j", as exposed in Section 5 and Appendix A. In Section 5.3, 
however, we deliver a variant of the existence proof of [16] for pure jump 
processes, which is possible in this case due to the absence of a diffusion 
component. 

Finally, Section 6 contains the proofs of the main results which build on 
the propositions of the previous sections. 

1.2. Notation. For the stochastic background and notation, we refer to 
standard text books such as [31] and [43]. We write = [0,oo) and M+_|- = 
(0, oo). Moreover: 

• Sfi denotes the space of symmetric d x d-matrices equipped with the scalar 
product {x,y) =Tr{xy). Note that Sd is isomorphic, but not isometric, to 
the standard Euclidean space We denote by < j} the 

standard basis of 5^, that is, the {kl)th component of c*-' is given by c^^ = 
Sik^ji + ^jk^iii^ — ^ij), where 5ij denotes the Kronecker delta. Additionally, 
we sometimes consider the following basis elements {e*-',i < j} which are 
positive semidefinite and form a basis of Sd- 



• St stands for the cone of symmetric d x d-positive semidefinite matrices, 
S'j"'' for its interior in Sd, the cone of strictly positive definite matrices. 
The boundary is denoted by dS'^ = S'^ \ Sj""*", the complement is denoted 
by {S'^Y, and S'^ U {A} denotes the one-point compactification. Recall 
that /Sj" is self-dual [w.r.t. the scalar product {x,y) = Tr(xy)], that is. 



Both cones, S'j' and S'J , induce a partial and strict order relation on Sd, 
respectively: we write x^yify — Sj", and x^yify — ^d~^ ■ 

• Md is the space of d x d-matrices and 0{d) the orthogonal group of di- 
mension d over JR. 

• Id denotes the d x d-identity matrix. 

Throughout this paper, a function / : S't; — t- R is understood as the restriction 
/ = dlsd of ^ function g : Md — )• M which satisfies g{x) = g{x~^) for all x E Md- 
Without loss of generality g{x) = f{{x + x~^)/2). We avoid using the vech 
operator, that is, to identify x £ Sd with a vector in M'^('^+^)/2 by stringing 
the columns of x together, while only taking the entries Xij with i< j. 

Throughout this article, we shall consider the following function spaces 
for measurable U C Sd- We write B{U) for the Borel c-algebra on U- bU 




if i=j, 
if i^j. 



S+ = {xeSd\{x,y)>0,yyeS+}- 
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corresponds to the Banach space of bounded real-valued Borel measurable 
functions f on U with norm ||/||oo = sup^g^ |/(rc)|. We write C{U) for the 
space of real- valued continuous functions / on U, Cb{U) for C{U) H bU, 
Cc{U) for the space of functions / G C{U) with compact support and Cq{U) 
for the Banach space of functions / G C(U) with lim^^-j^A f{x) = and norm 
ll/lloo = sup^g[/|/(x)| . Furthermore, C^{U) is the space of k times differen- 
tiable functions / on the interior of U, such that all partial derivatives of 
/ up to order k belong to C{U). As usual, we set C°°(f/) = HfeM C'^iU), and 
we write C^{U) = Cc{U) n C^{U) and C^{U) = Cb{U) n C'=(c7), for A: < oo. 

2. Definition and characterization of affine processes. We consider a 
time-homogeneous Markov process X with state space S'j' and semigroup 
{Pt)t>o acting on functions / E bS'^ , 

Ptfix) = [ f{OPt{x,dO, 

We note that X may not be conservative. Then there is a standard extension 
of the transition probabilities to the one-point compactification Sj" U {A} 
of Sj" by defining 

Pt{x, {A}) = 1 - pt{x, Sj), pt{A, {A}) = 1 

for all t and x G S'^ , with the convention that /(A) = for any function / 
on Sj". Thus X becomes conservative on Sj" U {A}. 

Definition 2.1. The Markov process X is called afftne if: 

(i) it is stochastically continuous, that is, liuis^tPsix, •) = pt{x, •) weakly 
on 5j~ for every t and x E Sj" , and 

(ii) its Laplace transform has exponential-affine dependence on the initial 
state 

(2.1) Pie~<"'^>=/ e~<"'«>pt(2;,(ie) = e-'^^*'"^-<^(*''')''^\ 

for all t and u,x £ S'j', for some functions (p '■ ^+ >^ S'^ ~^ ^+ and "0 '■ ^+ ^ 

Note that stochastic continuity of X implies that (j){t,u) and ijj{t,u) are 
jointly continuous in {t,u); see Lemma 3.2(iii) below. Moreover, due to the 
Markov property, this also means that pt{x,{A}) < 1 for all x € S'j" and 
t > 0. In contrast to [16], we take stochastic continuity as part of the defi- 
nition of affine processes, and consider the Laplace transform instead of the 
characteristic function. The latter is justified by the nonnegativity of X, the 
former is by convenience since, as we will see in Proposition 3.4 below, it 
automatically implies regularity in the following sense. 
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Definition 2.2. The affine process X is called regular if the derivatives 

dil){t,u) 



(2.2) F(n)^^'^(*'") 



dt 

exist and are continuous at n = 



R{u) 



t=o+ dt 



t=o+ 



We remark that there are simple examples of Markov processes which 
satisfy Definition 2.1(ii) but are not stochastically continuous; see [16], Re- 
mark 2.11. However, such processes are of limited interest for applications 
and will not be considered. 

In the following, we shall provide an equivalent characterization of the 
affine property in terms of the generator of X. As we shall see in (2.12), the 
diffusion, drift, jump and killing characteristics of X depend in an affine way 
on the underlying state. We denote hy x- some bounded continuous 

truncation function with xiO = C in 3' neighborhood of 0. Then the involved 
parameters are admissible in the following sense. 

Definition 2.3. An admissible parameter set {a,b, (3^^ ,c,j,m, fx) asso- 
ciated with X consists of: 

• a linear diffusion coefficient 

(2.3) a€S+, 

• a constant drift term 

(2.4) bh{d-l)a, 

• a constant killing rate term 

(2.5) ceR+, 

• a linear killing rate coefficient 
(2.6) 

• a constant jump term: a Borel measure m on S'j' \ {0} satisfying 



(2.7) / (llell Al)m(de)<oo, 

Js+\{0} 

• a linear jump coefficient: a d x d-matrix fj, = (pij) of finite signed measures 
on 5+ \ {0} such that fi{E) G 5+ for aU E G B{Sj \ {0}) and the kernel 

satisfies 

(2.9) / {x{Oiu)M{x,dC) <oo for aU x,^^^^ with = 0, 

Js+\{o} 
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a linear drift coefficient: a family (5^^ = (3^^ G Sd such that the linear map 
B : Sd ^ Sd the form 



(2.10) B{x) = Y,/3' 

satisfies 

(2.11) 



{B{x),u)- {x{0,^^)M{x,dO>0 
'sj\{0} 



for all X, n G with {x, u) = 0. 

We shall comment more on the admissibility conditions in Section 2.1 
below. The following three theorems contain the main results of this article. 
Their proofs are given in Section 6. First, we provide a characterization of 
affine processes on 5j~ in terms of the admissible parameter set introduced in 
Definition 2.3. As for the domain of the generator, we consider the space 5+ 
of rapidly decreasing C°°-functions on S'j", defined in (B.l) below. It is shown 
in Appendix B that e"<"'-> e S+, for u G 5++, as well as C^{Sj) C S+. 

Theorem 2.4. Suppose X is an affine process on S'^ . Then X is reg- 
ular and has the Feller property. Let A be its infinitesimal generator on 
Co(S'j'). Then 5+ C D{A) and there exists an admissible parameter set 
(a, b, (3"^^ , c, 7, m, fi) such that, for / G 5+ , 



d'fjx) 
2 Z-^ Qxij dxki 



+ B,,{x)) ^ - (c + (7, x))f{x) 



+ / {f{x + 0-f{x))m{dO 
'si\{o} 



(2.12) 



+ / (/(x + e)-/(x)-(x(0,V/(x)))M(x,dO, 
Js+\{o} 

where B{x) is defined by (2.10), M{x,d^) by (2.8) and 

(2.13) Aijki{x) = XikOji + Xiiajk + xjkau + xjiOik. 

Moreover, (j){t,u) andip{t,u) in (2.1) solve the generalized Riccati differential 
equations, for u G Sj", 

(2.14) ^^^ = F{i;{t,u)), HO,u)=0, 

(2.15) ^^!^ = R^^^t,u)), ^(0,n)=n, 
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with 

(2.16) F{u) = {b,u) + c- [ (e-<"'«> - l)?n((iO, 

Js+\{0} 



(2.17) 



R[u) = -2uau + B~^{u) + 7 

'e-<"'«>-l + (x(e),^) 



Sj\{0}\ llePAl 



;U((iO, 



w;/iere BJ-{u) = {f3^^ ,u). 

Conversely, let {a,b, /3^^ ,c,j,m, 1^1) be an admissible parameter set. Then 
there exists a unique afftne process on Sj" with infinitesimal generator (2.12) 
and (2.1) holds for all {t,u) € M4, x 5^", where (l){t,u) and il!{t^u) are given 
by (2.14) and (2.15). 

Remark 2.5. It can be proved as in [39] that X is conservative if and 
only if c = and ^p{t, 0) = is the only 5j'-valued local solution of (2.15) for 
u = 0. The latter condition clearly requires that 7 = 0. 

Hence, a sufficient condition for X to be conservative is c = and 7 = 
and 

/ miil^tjidC) + NjidC)) < 00 for all 1 < i < i < d, 

Js+n{\m>i} 

where fiij = fJ-fj — fJ-ij denotes the Jordan decomposition of fUj. Indeed, it 
can be shown similarly as in [16], Section 9, that the latter property implies 
Lipschitz continuity of R{u) on Sd. 

Due to the Feller property, as established in Theorem 2.4, any affine 
process X on S'^ admits a cadlag modification, still denoted by X (see, 
e.g., [43], Chapter III. 2). It can and will thus be realized on the space 
n = D(5+ U {A}) of cadlag paths a;:M+ ^ 5+ U {A} with w(s) = A for 
s > t whenever io{t—) = A or uj{t) = A. For every x € Sj", we denote by 
the law of X given Xq = x and by ( J-"^ ) the natural filtration generated by 
Xt. We also consider the usual augmentation 



(2.18) J^t:= n 



(x) 

t 



of (J-f^), where {J'j:^^) is the augmentation of (J't^) with respect to Fx- 

Then (J't) is right continuous and X is still a Markov process under (J't)- We 
shall now relate conservative affine processes to semimartingales, where^emi- 
martingales are understood with respect to the stochastic basis (^2, (J^)j , P^,) 
for every x. 
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Theorem 2.6. Let X be a conservative affine process on S"^ and let 
(a, 6, /3*-' , c = 0, 7 = 0, m, /i) be the related admissible parameter set associ- 
ated with the truncation function x- Then X is a semimartingale whose 
characteristics i/) with respect to x o-i"^ given by 

(2.19) At,ijki= [ Aijki{Xs)ds, 

Jo 

(2.20) Bt=f(b+[ x(.Om{dO + B{Xs))ds, 

Jo V ^5+\{0} / 

(2.21) u{dt, d^) = {m{di) + M{Xt,di)) dt, 

where B{x) is given by (2.10), Aijkiix) by (2.13) andM{x,d£,) by (2.8). Fur- 
thermore, there exists, possibly on an enlargement of the probability space, 
a d X d-matrix of standard Brownian motions W such that X admits the 
following representation: 

Xt = x + Bt + [ {y^sdWs S + E'^dWs y^) 
Jo 

(2.22) + f [ xmi^''{ds,di)-v{ds,dO) 

Jo Js+\{0} 

+ f I (e-x(e))/^''(d5,do, 

Jo Js+\{o} 

where S G satisfies S^S = a and denotes the random measure asso- 
ciated with the jumps of X . 

Hence, X is continuous if and only if m and /i vanish. 

Let V be the set of all families of probability measures (lPx)^g5+ on the 
canonical probability space {^,J-'^) such that {X , (¥" x) is a stochas- 
tically continuous Markov processes on 5^" with Px[^o = 2;] = 1, for all 
X G S'^ . Note that in contrast to [16], there is no need to impose regu- 
larity of X. For two probability measures P, Q on the convolution 
P * Q is defined as the push- forward of P x Q under the map (a;,a;') 1— )• 

Definition 2.7. An element (P^.) ^0+ G P is called: 

(i) infinitely decomposable, if for each k>l, there exists (^x^'^)x^s+ £ 'P 
such that 

p,.(i)+...+,.«=pS1)*---*pS); 
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(ii) infinitely divisible, if the one-dimensional marginal distributions Fx o 
X^^ are infinitely divisible, for all (t,x) G M+ x S'j'. 

In [16] it was shown that regular affine processes on x M" are infinitely 
decomposable Markov processes, and vice versa. In fact, this property was 
the core for the existence proof of affine processes in [16]. On S'^ the situation 
is different. The following counterexample reveals that not all affine processes 
on Sj" are infinitely divisible. 

Example 2.8. The affine process X on S'j' corresponding to the param- 
eter set (q = Id,b = 5Id, 0, 0, 0, 0, 0), where 6 £ [d — 1, oo), is the diffusion pro- 
cess initially studied by Bru [6]. By [6], Theorem 3, the Laplace-transforms 

E^[e-<^''">] = (det(/ + 2t'u))-^/2e-«^+2M-i«,x> 

are those of noncentral Wishart distributions WIS(5, d, x). By a well-known 
result due to Paul Levy, these Wishart distributions are not infinitely divis- 
ible if d> 2 (see [15], Section 2.C). 

Here, is our main result on infinite divisibility of affine processes on S'j". 

Theorem 2.9. Let d>2 and (IPx-)a;gs+ ^'P- The following assertions 
are equivalent: 

(i) (lPx)xg5+ infinitely decomposable. 

(ii) {X, (IF'z)2,.g5'+) ^"^ ^ffi''^^ with vanishing diffusion parameter a = 0. 

(iii) {X , {¥ x) is affine and infinitely divisible. 

2.1. Discussion of the parameters. We discuss and highlight some prop- 
erties of the admissible parameter set (a, 5, /S*-' , c, 7, m, /i) of an affine process. 
Let us therefore define the normal cone 

(2.23) Ng+{x) = {u£Sj\{u,x)=Q}, 

containing the inward pointing normal vectors, to S'j" at x € ."^ It will be 
shown in Lemma 4.1 below that Ng+{x) 7^ {0} only for boundary elements 

x^dSj. 



^Indeed, we obtain (2.23) from the general definition in (A. 8) below by choosing y = 
and y = 2x, and using the self-duality of : {u, y) >0 for all y,u£ . 
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2.1.1. Diffusion. The diffusion term does not admit a constant part, and 
its linear part is of the very specific form 

{u, A{x)u) = 4(x, uau). 

This property of A{x) has also been stated in the setting of symmetric cones 
in [26]. We could thus write the second order differential operator in (2.12) 
as 

5E-4««W^ = 2(..VaV/(.)). 

The reason why we introduce and use the symmetrization (2.13) of A{x) is 
that it corresponds to the quadratic characteristic (2.19) of the semimartin- 
gale X. 

2.1.2. Drift. The remarkable drift condition (2.4) has been assumed in 
many previous papers. Here is the first time where necessity and sufficiency 
of (2.4) are proved in the full generality in the presence of jumps. Note 
that in dimension d=l, the drift condition simply reduces to nonnegativity 
6 > 0. But for dimension d>2, the boundary of the state space S'^ becomes 
curved and kinked, implying a nontrivial trade-off between diffusion a and 
b. 

Concerning the form of B, let us note the following: condition (2.11) 
implies in particular 

(2.24) /3(^^,^ - / ^|^^,,(de) G Si, for ah 1 < . < d, 

Js+\{o} ll?lr ^ ^ 

where for any matrix u £ Sd, u\{i} denotes the matrix where the ith row 
and column are deleted. Indeed, inserting x = c" in condition (2.11) yields 

{B{cn,u)- [ j^|l^^,,(dO>0 

Js+\{o} ll?lr ^ ^ 

for all u £ S'j' with {d'^^u) = 0. Since the ith column and row of such an 
element u G Sj" is zero, it follows that 

{B{c%u)- [ j^|)i^^,,(de) 

Js+\{o} IkII Al 

(2.25) 

- - 4,^^^^ iiePAi ^^^^'^^ ^ 

By choosing appropriate elements u\{i} G '^d-i^ '^^ can further derive the 
integrability of xiOki for all k ^ i, which implies 

Js+\{0} a Al \j5+\{0} C Al 
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As (2.25) and (2.26) must hold true for all u\{i} G 'S'JLi, assertion (2.24) is 
proved. 

Note that the (ij)th component of the adjoint operator is given by 

(2.27) Bj^iu) = {/3^^,u), 

since {B{x),u) = 1^'^ Xij,u) = J2^j{/3'^ ,u)xij = {B'^{u),x). 
In most previous papers, B{x) is of the form 

(2.28) B{x) = Hx + xH~^. 
In this case, 

(2.29) {B{x),u) = {Hx + xH'^ ,u) =0 for all n G 5j" with (x, n) = 0, 
and hence (2.11) is equivalent to 

/ {x{0,^)M{x,dO=0, 
Js+\{0} 

for all X, u G 5^" with {x, u) = 0. 
If B{x) is of the form 

(2.30) B{x) = Hx + xH'^ +T{x), 

where H G Md and F : 5^ — )■ 5^ linear satisfying r(5j') C then, in view 
of (2.29), condition (2.11) holds true as long as 

(r(x),n)-/ {x{i),u)M{x,di)>Q 

for all x,?x G with {x,u) = 0. As a bold conjecture, we claim that any 
B{x) satisfying (2.11) is of form (2.30). 

Here is a simple example where B{x) is of the form (2.30) but not of the 
usual form (2.28): let d = 2 and 

\Xi2 Xii J 

It can be easily checked that (2.11) is satisfied, while B{x) cannot be brought 
into the form (2.28). If xu models the (squared) volatility of the iih. stock 
price, as in (1.1), then this drift specification admits level impacts of the 
volatility of stock 1 on the volatility of stock 2, and vice versa. 



2.1.3. Killing. See Remark 2.5. 
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2.1.4. Jumps. Condition (2.7) means that jumps described by m, which 
can for instance appear at x = 0, should be of finite variation entering the 
cone Sj", since infinite variation transversal to the boundary would let the 
process leave the state space. Similarly, condition (2.9) asserts finite variation 
for the inward pointing directions, while we could a priori have a general 
jump behavior (supported by Sj" due to the affine structure) parallel to 
the boundary. Note that in the case d=l, which corresponds to R+, the 
linear jump part can have infinite total variation (see [16], equation (2.11)). 
However, due to the geometry of the cone S'j", we conjecture that in higher 
dimensions d>2 such a behavior is no longer possible and that all jumps are 
in fact of finite total variation. In any case, for d>2, affine positive matrix 
valued diffusion processes cannot be approximated (in law) by pure jump 
processes, since this would yield a contradiction to condition (2.4). See also 
Remark 5.12 below. 

3. Affine processes are regular and Feller. Suppose X is an affine process 
on Sj". The main result of this section is that X is regular in the sense of 
Definition 2.1. In addition, we shall prove that Pt is a Feller semigroup on 
Co(S'j'). In order to show both properties, we shall mainly rely on Lemma 3.3 
below. The Feller property is then a simple consequence of this statement 
and regularity is obtained by arguing as in Keller-Ressel, Schachermayer 
and Teichmann [35], who obtained the corresponding statements for affine 
processes on the state space x M"; see [35], Theorem 4.3, and also the 
Ph.D. thesis of Keller-Ressel [34]. We observe that most arguments of [35] 
translate to our setting without major changes. It is only required to tailor 
some technicalities to the cone 5^". We start with the following elementary 
observations. 

Lemma 3.1. Ifu£ dS^ and S'^ Bv ^u, then v G dSj. 

Proof. Let x G S'j" \ {0} such that (x, n) = 0. Then, Bv implies 
0<{v,x)<{u,x) = O.Rence,v£dSj. □ 

We now derive some first properties of the functions (j) and ip in (2.1). 

Lemma 3.2. Let X be an affine process on Sj". Then, we have: 
(i) The functions (j) and ip satisfy 



(3.1) 
(3.2) 



(j){t + s,u) = (j){t, u) + (j){s, ■tp{t, u)) 
'4}{t + s,u) = ip{s,ip{t,u)) 



for all t, s G E+. 
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(ii) For all u,v € Sj" with v and for all t>0, the order relations 
(3.3) 4>{t,v) < (l){t,u) and ijj{t,v) ^ ilj{t,u) 

hold true. 

(iii) The functions 4> and are jointly continuous in 1R+ x Sj" . Further- 
more, (l){t,u) and u>-^ tpitju) are analytic on S'^'^ . 

Proof. Assertion (i) follows directly from the Chapman-Kolmogorov 
equation, 

Js+ Js+ 

^ ^-<P{t,u)-<l>{s,-4,it,u))-{^P{s,^P{t,u)),x) _ 

For the proof of (ii), note that v is equivalent to {v,x) < {u,x) for all 
X G Sf. By the monotonicity of the exponential function, we have for all 
X E and for all t>0, 

^-4>it,v)-m,v),x) ^ I e"^''^^^pt{x,dC)> [ e-<"'«>pt(2;,(ie) 
■^s+ Js+ 

— g-<l>{t,u)-{i){t,u),x) 

and the assertion follows by taking logarithms. 

Concerning statement (iii), note that stochastic continuity of X implies 
joint continuity of P^e"^"'^^ in {t,u) G M-|_ x S'^ (this follows, e.g., from [4], 
Lemma 23.7), for all x G S'j'. This in turn yields continuity of the func- 
tions {t,u) I— )■ <j){t,u) and {t,u) i— )• 'il){t,u). The second assertion follows from 
analyticity properties of the Laplace transform. □ 

The following property of 'ijj is crucial. 

Lemma 3.3. Let il: -M^ x Sj ^ Sj he any map satisfying ip^Oju) = u 
and the properties (i)-(iii) of Lemma 3.2 (regarding the function tp). Then 
ij{t,u) e for all {t,u) G M+ X S'j"+. 

Proof. We adapt the proof of [34], Proposition 1.10, to our setting. 
Assume by contradiction that there exists some (t, u) G M+ X 5^+ such that 
ip{t,u) £ dS'^ . Let us consider the interval (0, Xmmiu)] / 0, where Xmin{u) > 
denotes the smallest eigenvalue of u. Then for all v G (0, Amin(^*)] we have 
vid ^ u- Since 'ip{t,u) admits property (ii) of Lemma 3.2, we obtain 

S+^^p{t,vId)^Ht,u)e^S+ 
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for all V G (0, Amin(^i)]- Consequently, Lemma 3.1 yields that ilj{t,vld) G 
Hence, 

det{^P{t,vId)) = 

for all V G (0, Amin(^)] • The analyticity of n i— tlj{t, u) on S'j""'" carries over to 
u I—)- det(V'(i, It)) and implies that det(^/'(t, v/^)) = for all v G Indeed, 
the set of zeros of det{tjj{t,vl(i)) has an accumulation point in M+_|_, which 
implies that det(V'(t, t^/rf)) vanishes entirely on ]R+_|_. The same statement 
holds true for t replaced by |. Indeed, if S f^'S'j'; then the assertion 

is shown by the same arguments as above. Otherwise, if iIj{^,u) G we 
have for all v G with vId ^ that is, for all v G (0, Amin('0(f i ^)] 

which yields again ilj{^,vld) G OSj" and det(?/'(|, f I^)) = for all v G (0, 
Amin('0(|) ^)]- The same reasoning as before then leads to det(V'(|, f/d)) = 
for all V G By reapplying this argument, we finally get for every n G N 
and for all v G 

From the continuity of the function 1 1— t- ■ilj{t,u) and of the determinant, we 
deduce that for any v G , 

= Jim det (^^ (^-^,vld^ ) = det(V(0, vId)) = det{vld) = v'^>0, 

a contradiction, and the assertion is proved. □ 

We may now formulate the main result of this section. 

Proposition 3.4. Let X be an affine process with state space Sj". Then, 
we have: 

(i) X is a Feller process. 

(ii) X is regular. 

Proof. In order to prove (i), it suffices to show that for all / G Co(S'j') 

(3.4) ^hm^Pt/(x) = fix) for aU x G S+ , 

(3.5) Ptf G CoiSj) for ah t G M+, 
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(see, e.g., [43], Propostion III. 2. 4). Property (3.4) is a consequence of stochas- 
tic continuity, which imphes for ah / G Co{S'^) and x G Sj" 

hm Ptf{x) = f{x). 

Concerning (3.5), it suffices to verify this property for a dense subset of 
Co(<S'j~). By a locahy compact version of Stone- Weierstr ass' theorem (see, 
e.g., [48]), the linear span of the set {e-<"'^> | u G is dense in Co{Sj). 

Indeed, it is a subalgebra of Co(<S'j'), separates points and vanishes nowhere, 
as all elements are strictly positive functions on Sj". From Lemma 3.3, we 
can deduce that Pje"^"'^^ G Co{Sj) if u G Sj""^, since '4'it,u) G and 
{^p{t, u), x) > for X / implying that 

p^g-(«.2:> — Q-<l>it,u)-{tp{t,u),x) 

goes to as x —7- A. Hence, statement (i) is proved. 

The proof of (ii) follows precisely the lines of [35], proof of Theorem 
4.3. Using Lemma 3.3, one may mimic the proof of [35], Theorem 4.3, to 
obtain that differentiability of ■0(t, u) inu £ •S'j'^, which follows from Lemma 
3.2(iii), implies differentiability of 'il){t,u) in t for t = and for all u G S'^ . 
□ 



By the regularity of X, we are now allowed to differentiate the equations 
(3.1) and (3.2) with respect to t and evaluate them at t = 0. As a conse- 
quence, (j) and satisfy the system of differential equations 

^'^(*'^)^F(V^(t,^)), ,^(0,n)=0. 



dt 

d'4){t,u 



Rm,u)), V(0,^x)=nG5+, 



where F and R are defined as in (2.2). The analysis of these (generalized 
Riccati) differential equations is subject of Section 5.1, whereas the specific 
form of F and R is elaborated in the following. 



4. Necessary parameter restrictions. In this section, we derive necessary 
parametric restrictions, that is, given an affine process on S'j', we determine 
necessary implications on a set of parameters which only ensue from Defi- 
nition 2.1. These conditions are precisely the conditions on the admissible 
parameter set as of Definition 2.3. The form of the functions F and R as 
defined by (2.2) is then characterized by means of this parameter set, which 
is stated in Proposition 4.9 below. For its proof, we first provide a number 
of technical prerequisites. 
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Lemma 4.1. Letx,u£S^and 

(4.1) x = OAO^ = Odiag(Ai >0,...,Ad_r >0,0,...,0)O^ 

be the diagonalization of x with r > and O € 0{d). Then the following 
assertions are equivalent: 

(i) ux = xu = 0, 

(ii) {x,u) =0, 

(iii) u is of form 

(4.2) u = 0(l l)0- 



with w G St . 



Proof. The direction (i) =^ (ii) is obvious. In order to prove the imph- 
cation (ii) =^ (iii), define v as v = O^uO. Then we have 

{)={x,u) = {k,6^u0)= ^ \iV^^, 

i<d—r 

which imphes vu = for alH < d — r and by the positive definiteness of v it 
must then be of form 


w 

with w G . Thus, u is given by (4.2). This then imphes that ux = xu = 0, 
which proves the direction (iii) =^ (i). □ 

Lemma 4.2. Letp be an orthogonal projector, that is, p G S'j" and p"^ =p 
(see, e.g., Kato [33], Section I.6.7J, and define q = Id—p. Then q is an 
orthogonal projector and the orthogonal complement of p in S^ equals 

{veSj\ {p, v)=0} = {quq \ u G Sj}. 

Proof. That q is an orthogonal projector fohows by inspection. The 
diagonalization of p is of the form p = OAO^ with A = diag (1,...,1,0,...,0), 
and thus q = 0{ld - A.)0'^ . In view of Lemma 4.1, we conclude that v G S'j' 
is orthogonal to p if and only if v = qvq. This proves the assertion. □ 

Lemma 4.3. Let u be in Sd and x G dS'^ such that ux = xu = 0. Then, 
the linear map defined by 

Tu'.Sd^ Sd, V ^ TuV := uvu 

has the following properties: 
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(i) Tu is self-adjoint and Tu{S'^ — R+x) C Sj". 

(ii) There exists an element v £ Sd such that T^v = u. 

Proof. The assertion (i) is obvious, since for every k G M_|_, T^kx = 
kuxu = and T^v = uvu G S'j" ii v £ S'^ . For proving part (ii), we use the 
fact that X is of form (4.1) and that all zero divisors u in of x can be 
represented by (4.2) with w € Sr- Thus, setting 

where satisfies ww~^w = w, yields T^v = u. □ 

Lemma 4.4. Let V denote a vector space^ over M. Let L:S'^ -^V he 
an additive (resp., homogeneous additive) map, that is, for all x,y G Sj" and 
A = 1 (resp., for all A G M+ ) we have 

L{x + Ay) = L(x) + \L{y). 

Then L{x) is the restriction of an additive (resp., M.-linear) map on Sd. 

Proof. We define the map L:Sd^V as 

L{x -y) := L{x) - L{y), x,yeSj. 

L is well defined, as for u, v,x,y ^ S'j' such that u — v = x — y we have 

L(n) - L{v) = L{u -v) = L{x - y) = L{x) - L{y). 

Since S'j' — = Sd, the domain of L is all of Sd. Also, L(0) = by the 
additivity of L. Hence, L is the restriction of L to Sj". Homogeneity of L 
holds, as for A > 0, z = X — y G 5(i we have by definition 

L{\z) = L{Xx) - L{\y) = XL{x) - \L{y) = \L{z). 

Finally, we show additivity of L. Choose w,z £ Sd such that^z = x — y,w = 
u — V, hence w + z = (x + u) — {y + v). By the definition of L, we have 

L{z) = L{x)-L{y), L{w)=L{u)-L{v), 

and by the additivity of L we obtain 

L{w + z) = L{x + u)- L{y + v) = L{x) + L{u) - L{y) - L{v) 

= L{z) + L{w). □ 



*In the proof of Proposition 4.9 below, V corresponds to Sd, the vector space of hnear 
maps Sd, or the vector space of finite signed measures on Sd- 
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We now provide a convergence result for Laplace transforms (in fact 
Laplace-Fourier transforms) , which is most relevant for the analysis of affine 
processes. 

Lemma 4.5. Let f„ he a sequence of measures on Sd with 

Ln{u) = / e'^^'^^t'nfd^) < oo and lim Ln(u) = L(u), uGSt, 
Jsa 

pointwise, for some finite function L on S'^ continuous at u = 0. Then Un 
converges weakly to some finite measure v on Sd and the Fourier-Laplace 
transforms converge for u G S'^'^ U {0} and v £ Sd to the Fourier-Laplace 
transforms of v, that is, 

lim [ e-<"+^^'«>i/„(dO= f e-^^'+'^'^^uidO- 
Jsd Jsd 

Ln particular, u{Sd) = lim„_>.oo Vn{Sd) and 



L{u)= / e-<"'«>i/((ie), 



for a//uG 5j"+U{0}. 



Remark 4.6. Instead of n = we could take any set K of points at the 
boundary K C dS'^: if we assume continuity of L at points in K, then we 
obtain the equality of L with the Laplace transform of u for all points in 
K. Additionally, continuity is too strong an assumption, since we only need 
right continuity of L along the segment u + eld for e = at the points from 
the boundary under consideration. 

Proof of Lemma 4.5. Since Vn{Sd) = Ln{Q) is bounded, we know by 
general theory that Vn has a vague accumulation point v, which is a finite 
measure on Sd- 

Since Ln{u) < cxd on Sj", it follows by well-known regularity properties of 
Laplace transforms (see, e.g., [19], Lemma 10.8) that the functions L^ admit 
an analytic extension on the strip Sj"*" + iSd, still denoted by L„: 

{u + iv)^Ln{u + iv)= I e-<"+™'«>i/„(dO. 



Moreover, pointwise convergence of the finite convex functions L„ to L on S^ 
implies that this convergence is in fact uniform on compact subsets of 5j" 
(see, e.g., Rockafellar [44], Theorem 10.8). Hence, the functions L„ are uni- 
formly bounded on compact subsets of S^'^ and since \Ln{u + iv)\ < Ln{u), 
also on compact subsets of Sj"*" -|- iSd- Therefore, and since Sj""*" is a set 
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of uniqueness in S'j""*" + iSd, it follows by Vitali's theorem ([41], Chapter 1, 
Proposition 7) that the analytic functions L„ converge uniformly on com- 
pact subsets of S'j""'' + iSd to an analytic limit thereon. By Levy's continu- 
ity theorem, we therefore know that for any u € S'^^ the finite measures 
exp(— (u, converge weakly to a limit, which by uniqueness of the 
weak limit has to equal exp(— Whence the only vague accumu- 
lation point of f„ is i^. Vague convergence implies weak convergence if mass 
is conserved. Continuity of L{u) at ti = implies this mass conservation: 
indeed, by weak convergence of e~^'^^'*'^^i'„ we arrive at 

L{eld)=lim [ e-<"^'''«V„(dC) = / e-^'^'"^^ iy{dO 

and therefore — by dominated convergence — we obtain that the limit e — )• 
yields 

L(0)= / u{dO, 
J Si 

which is the desired mass conservation, hence weak convergence, which 
means in turn convergence of the Fourier-Laplace transform at ti = 0. □ 

Finally, let us state a general comparison result for ODEs and hereto 
introduce the notion of quasi-monotonicity, which we shall need several times 
throughout this article, in particular in the proofs of Propositions 4.9 and 
5.3 below. 

Definition 4.7. Let U C Sd he an open set. A function f:U^Sd is 
called quasi-monotone increasing if for all elements x,y £U, n € S'J' which 
satisfy x^y and {x,u) = {y,u), 

{fix),u)<{fiy),u) 

holds true. Accordingly, we call / quasi- constant if both / and — / are quasi- 
monotone increasing. 

The following comparison result can be deduced from a more general 
theorem proved by Volkmann [52]. 

Theorem 4.8. Let U C Sd be an open set. Let f -.[OjT) x U ^ Sd be a 
continuous locally Lipschitz map such that f{t,-) is quasi-monotone increas- 
ing on U for all t £ [0, T). Let <tQ <T and x,y: [0,to) U be differen- 
tiable maps such that x{0) ^ y(0) and 

x{t) - f{t, x{t)) ^ y{t) - fit, y{t)), 0<t<to. 

Then we have x{t) ^ y{t) for all t G [0,to)- 
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4.1. The functions F and R. The main result of this section characterizes 
the form of the functions F and R as defined by (2.2). 

Proposition 4.9. Let X be an afftne process with state space Sj". Then 
there exist parameters (a,6,/3*-^,c,7,m,/i), where a, 13^^ ,c,j,m, fj, satisfy the 
admissibility conditions of Definition 2.3 and h G S'^ , such that the functions 
F and R are of the form (2.16) and (2.17). 

Remark 4.10. Note that for the moment we only obtain b G Sj", and 
not (2.4). 

Proof of Proposition 4.9. As the proof of Proposition 4.9 is rather 
long, we divide it into several steps: 

Step 1. Necessary admissibility conditions for b,c,^,m. In order to derive 
the particular form of F and R with the above parameter restrictions, we 
follow the approach of Keller- Ressel [34], Theorem 2.6. Note that the t- 
derivative of P^e"^"'^^ at t = exists for all x,u £ , since 

p^Q-{u,x) _ ^-{u,x) ^-4>{t,u)-{-tp{t,u),x) _ ^-{u,x) 

lim = lim 

t^o+ t t^0+ t 

(4.3) 

= (-F(tx)-(i?(n),x))e-<"'^> 
is well defined by Proposition 3.4. Moreover, we can also write 

-F{u) - {R{u),x) 

Pte'^"^'^^ — e~^"'^^ 
= lim — ^ 

= lim - 

t->o+ t 



->o+ t \Js+\{0} J 
lim f ^ / (e-<-'«) - l)p,{x, dC + x) + PJ^^l^ll^) . 



By the above equalities and the fact that pt{x,S'^) < 1, we then obtain for 
^ = 

0> lim ^iM£l^ = _i7(o)-(i?(o),x). 

t-^o+ t 

Setting F{0) = c and -R(O) = 7 yields c E as in (2.5) and 7 G Sj" as in 
(2.6). We thus obtain 

-iF{u)-c)-{Riu)--f,x) 

(4.4) 

= lim - /" (e-<"'«> - l)pt(2;,de + a;). 
t^o+ t Js+-x 
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For every fixed t > 0, the right-hand side of (4.4) is the logarithm of the 
Laplace transform of a compound Poisson distribution supported on 5'j' — 
R+x with intensity pt(x,5j~)/t and compounding distribution pt{x,d^ + 
x)/pt{x, S'^). Concerning the support, note that the compounding distri- 
bution is concentrated on S'^ — x, which implies that the compound Poisson 
distribution has support on the convex cone Sj" — M+x. By Lemma 4.5, the 
pointwise convergence of (4.4) for t — )• to some function being continu- 
ous at 0, implies weak convergence of the compound Poisson distributions 
to some infinitely divisible probability distribution K{x, dy) supported on 
S'^ — M+x. Indeed, this follows from the fact that any compound Poisson 
distribution is infinitely divisible and the class of infinitely divisible dis- 
tributions is closed under weak convergence ([47], Lemma 7.8). Again, by 
Lemma 4.5 the Laplace transform of K{x, dy) is then given as exponential 
of the left-hand side of (4.4). 

In particular, for x = 0, i^(0, dy) is an infinitely divisible distribution with 
support on the cone . By the Levy-Khintchine formula on proper cones 
(see [49], Theorem 3.21), its Laplace transform is therefore of the form 

expf-(6,n) + / (e-<"'«> - l)m{di)\, 
V Js+\{0} J 

where h G Sj" and m is a Borel measure supported on S'^ such that 

/ (lieil Al)m(fie)<oo, 

yielding (2.7). Therefore, 

F{u) = {b,u)+c- [ (e-<"'«> - l)m{dO. 
Js+\{0} 

Step 2. Necessary admissibility conditions for /3*-',/x. We next obtain the 
particular form of R. Observe that for each x G S'^ and /c E N, 

exp(-(F(n) - c)/k - {R{u) - 7, x)) 

is the Laplace transform of the infinitely divisible distribution K{kx, dy)*^^^ , 
where *^ denotes the ^ convolution power. For A; — )• 00, these Laplace 
transforms obviously converge to exp(— (i?(u) — 7,x)) pointwise in u. Us- 
ing again the same arguments as before [an application of Lemma 4.5 as 
below (4.4)], we can deduce that K{kx^dy)*^^^ converges weakly to some 
infinitely divisible distribution L{x,dy) on S'^ — R+x with Laplace trans- 
form exp(— (i?(u) — 7, x)) for u G Sj". 

By the Levy-Khintchine formula on 5*^ ([47], Theorem 8.1, indeed on 
]^{d(d+i)/2) -^y jxiodifying the scalar product appropriately), the characteristic 
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function of L{x,dy) has the form 

L{x, u) = exp ^-(n, A{x)u) + {B{x),u) 



Sd\{0} 



for u G iSd, where A{x) is a symmetric positive semidefinite hnear operator 
on Sd, B{x) G Sd, M{x, •) a measure on Sd \ {0} satisfying 



A l)M{x,dO < oo, 

Sd\{0} 

and X some appropriate truncation function. Furthermore, by [47], Theorem 
8.7, 



Is„ 



f{0lpt{x,d£,+x 
(4.5) ' 



f{Om{dO+ [ f{OM{x,dO 
Sd\{0} JSd\{0} 



holds true for all f : Sd^^ which are bounded, continuous and vanishing 
on a neighborhood of 0. We conclude that M{x,d^) has support in — x. 
Therefore, the characteristic function L{x,u) admits an analytic extension 
to S'^ X iSd, which then has to coincide with the Laplace transform for 
u S . We conclude that, for all x G Sj", 

-{R{u) -7,x) 
(4.6) =^{u,A{x)u)-{B{x),u) 

+ I (e-<"'«>-l + (x(0,"»M(x,dO, uGSj. 

As the left-hand side of (4.6) is linear in the components of x, it follows 
that X H> A{x), x H> B{x) as weh as x i-)- /gjdlClP A l)M{x,d(,) for every 
E G B{Sd \ {0}) are homogeneous additive maps on S'j' in the sense of Lemma 
4.4. This then implies that they are restrictions of linear maps on 5^, such 
that we can write 



A{x) = '^aijXij, B{x) = 
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where (recall that &^ denotes the standard basis of Sd defined in Section 
1.2): 

aij = Qji = (1 + 6ij) — - — :Sd^Sd linear, 

and 

E ^ Hij{E) = fiji{E) = (1 + Oij)— 

are finite signed measures on S'^^ \ {0}. The fact that M{x,-) is a nonneg- 
ative measure for each x E Sj" implies immediately that fJ.{E) is a positive 
semidefinite matrix. 

In (4.5), take now x = ^e*-' and nonnegative functions f = fn ^ Cb{Sd) 
with /„ = on Sj" — ^e^-' ■ Then for each n the left-hand side of (4.5) is zero 
since the Pt(^e*-', + ^e*-') is concentrated on — ;^e*-'. As supp(?n) C 5j~, 
the first integral on the right vanishes as well. Hence, 

0= / m)M(-e^\dc) = [ ^^/-e^^/.(dO 

= - / mrl^(^'^i(dO + (1 - s,J){^^,,{dC) + 2fiij{dm 

^ JSa\{0} ll?lr ^ ^ 

for any nonnegative function /„ G Ch{Sd) with = on 5j~ — ^e*-^ implies 
that supp(^,j) C Sj" — ^e"^^ for each n. Thus, we can conclude that supp fiij C 
5+ for ah l<i,j<d. 

Now let T : Sd ^ Sd be any linear map with the property T^S'^ — M+x) C 
Sj". Then T(supp(L(x, dy))) C 5^. This implies that the pushforward T^:L{x, •) 
of L{x,dy) under T is an infinitely divisible distribution supported on S'j". 
By the Levy-Khintchine formula on proper cones (see [49], Theorem 3.21, 
and by [47], Proposition 11.10) this implies that for all x G S'j' 

(4.7) TA{x)T^ = 0, 

(4.8) TB{x) + / (x(re) - TixiOmix, d^ S 5+, 

Js+\{0} 



(4.9) / (llell Al)r,M(x,de)<oo, 

Js+\{0} 

where x denotes some truncation function associated with T^L{x,-) and 
r^,M the pushforward of M under T. Due to (4.9), we can set x = 0. Thus, 
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(4.8) becomes 

(4.10) TB{x) - I T{x{i))M{x, di) G 5+. 

Moreover, equations (4.7), (4.10) and (4.9) are equivalent to 

(T^u, A{x)T'^v) = for all v G 5^, 

{B{x),T^v)- [ {{x{0),T^v)M{x,dO>0 forallt;G5+, 
Js+\{0} 

(4.11) / (llrell Al)M(x,dO<oo. 



In particular, we claim that 

{u,A{x)u) = 

(4.12) 

for all u £ Sd s.t. ux = xu = 0, 



(4.13) 

for all u G S'^ s.t. ux = xu = 0, 



{x{0,u)M{x,dO<oo 

(4.14) 



for all u G s.t. ux = xu = 0. 

Indeed, if x is invertible then urr = is equivalent to u = and the assertions 
are obvious. Otherwise, if x is in dS'^ , the linear map T„ defined in Lemma 
4.3 is self-adjoint and satisfies Tu{S^ — M+x) C Furthermore, by Lemma 
4.3(ii), there exists an element v £ Sd such that T^v = u. Hence, 

{u,A{x)u) = {tJv,A{x)tJv)=0. 

It follows from the proof of Lemma 4.3 that for u G S*^ , ^; is an element of Sj" 
as well and we have {B{x),u) = {B{x),TJ v) and (x(0''") = {{x{C))iT~^ v). 
Equation (4.14) is obtained by choosing T = T^ in (4.11). Indeed, 



{^,u)M{x,dO= / {Id,Cu)M{x,dO 
5+n{||Cil<i} ^s+n{i|5i|<i} 



<||/,|| / Uu\\M{x,dO 
■Js+n{m\<i} 



\\T^aM{x,dO <oo. 
s+n{||5||<i} 
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From these arguments and Lemma 4.1, properties (2.11) and (2.9) can be 
derived so far. Thus, only (2.3) remains to be shown. 

Step 3. Necessary admissibility condition for a. Due to the hnearity of 
A{x), {u,A{x)u) can be written as 4(x,'!?(n)), where the (ij)th component 
of '&{u) E Sd is defined by "Qij^u) = 1/A{u,aiju). Note that -& is defined on 
ah of Sd- Given that for all x G S'^ , A{x) is a positive semidefinite operator 
on Sd, {u,A{x)u) > for all Sd and therefore, by the self duality of S'j", 
'd{u) G Sj". By (4.12), we have for all u such that ux = xu = 

(4.15) 0= {u,A{x)u) =4{x,'d{u)). 

Next, we show that i9 is quasi-constant, that is, {x, t!){u + w) — i?(ti)) = for 
all x,u,w ^ Sj" with (x, w) =0 (see Definition 4.7). Indeed, pick x,u,w ^ 5j~ 
with {x,w) = 0. According to our assumptions, A{x)w = 0, due to (4.15) and 
the positivity of A. Hence, 

4(x, '&{u + w) — 'i?(ii)) = {u + w, A{x){u + w)) — {u, A{x)u) 

= {u, A{x)w) + {A{x)w, u) = 0, 

where the second last equality holds in view of the symmetry of A{x). 

We now claim that there exists some a S Sj" such that •d{u) = uau, for 
each u £ Sd- It is sufficient to show that this statements holds for all or- 
thogonal projectors p € Sj", that is, there exists some a € S'j' such that 
'&{p) =pap for all orthogonal projectors p. Indeed, if this is the case, we can 
derive the general statement in the following way: take u G Sj", then — by 
spectral decomposition — there are numbers Aj > and orthogonal projec- 
tors Pi, which are mutually orthogonal, such that u = AjPi (see, e.g., 
Kato [33], Section 1.6.9). Since the assertion holds for all orthogonal projec- 
tors, we have that 

d 

2^{u) = J2 KXMm+Pj) - ^Pi) - ^{Pj)) 

by the property that is quadratic. Since pi + pj is again an orthogonal 
projector, we obtain the result. 

We prove the assertion on orthogonal projectors by quasi-constancy. Take 
an arbitrary orthogonal projector p and define q = Id — P- Additionally, we 
define a = 'd{Id)- By quasi-constancy, we obtain 

(x, ^{p + q)- ^{q)) = {y, {){p + q) - d{p)) = 

and 

(x,T9(p)) = (y,T?(g))=0, 

for all x,y G S*^ with (x,p) = and {y,q) = 0. Therefore, a — ^{q) and §{p) 
are orthogonal to the orthogonal complement oi p in S^ (i-e., the positive 
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symmetric matrices of the form quq by Lemma 4.2), and a — 'd{p) and '&{q) 
are orthogonal to the orthogonal complement of q in Sj" (the positive sym- 
metric matrices of the form pup by Lemma 4.2). This means that we can 
write 

a = {){p)+{){q)+P, 

where the symmetric matrix j3 is orthogonal to all elements which are or- 
thogonal to p and q (in S^), that is, /3 is orthogonal to the linear span of 
matrices of the form pup and quq. However, such a decomposition is unique, 
since all vectors in the sum are mutually orthogonal, and the decomposition 
is given by 

a= {p + q)a{p + q) = pap + qaq + {paq + qap). 

Therefore, we can conclude the assertion i!){p) =pap. Since p was arbitrary 
the assertion is proved. 

Finally, all the derived restrictions on the parameters together with (4.6) 
then yield (2.17). □ 

Remark 4.11. An alternative proof for the special form of the diffusion 
matrix A{x) can also be established by Stokes' theorem [50] on Riccati 
ODEs. 



4.2. Infinitesimal generator. The aim of this section is to prove the form 
of the infinitesimal generator as stated in (2.12). 

Proposition 4.12. The infinitesimal generator A of an affine process 
on S'j' satisfies 5+ C T^iA) and is of the form (2.12) for all f G 5+ and 

Proof. As already mentioned in the proof of Proposition 4.9, the t- 
derivative of P^e"^"'^^ at t = exists pointwise for all x,u£ Sj" and is given 
by (4.3). Furthermore, x ^ (-F(u) - (i?('u), 2;))e-<"'^> G Co{Sj), for u G 
S'^^ . As (Pt) is a Feller semigroup on Co(5'j'), it follows from [47], Lemma 
31.7, that {e-<"'^> | u £ 5++} G D{A) and 

_^g-(n,x> ^ (^-F{u) - (i2(n),x))e-<"'^>. 

Combined with Proposition 4.9, we thus obtain 

^e-<"'^> = (-{b,u)-c+ [ (e-<"'«> - l)m(dO 
V ^s+\{0} 

+ lluau- B^{u) -7 
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'e-<"'«)-l + (x(e),^) 



+ 



SI\{0} 



iiepAi 



(4.16) 



ijki[x)UijUkie + {b + B{x),Ve 



i,j,k,l 

-(c+(7,x))e-<"'-) 



+ 



+ 



S,+\{0} 



s2\{o} 



Indeed, in order to obtain the form of the diffusion part, observe that we 
have by symmetrization 



2{uau,x) = 2 ^ ajkXiiUijUki 



i,j,k,l 



^ ^ {xikOiji + Xiiajk + XjkOtii + Xjiaik)uijUki 

i,j,k,l 

\ Ajkl{x)uijUM; 
i,j,k,l 



see (2.13). 



According to Theorem B.3, the linear hull A4 of {e 



-in,-) 



U 



is 



dense in 5+ with respect to the family of seminorms Pk,+ defined in (B.2). 
Denoting the right-hand side of (2.12) by A", we now claim that for every 



(4.17) 



hm P»/n-^»/||oo=0, 



where /„ G 7W such that lim„_j.ooPfc,+(/ — fn) = for every k. Indeed, this is 
obvious for the differential operator part of A". By choosing x{€) = l{||C||<i}C 
and by denoting g{x) := fn{x) — f{x), we obtain the following estimate for 
the integral part: 



s:\{o} 



g{x + - g{x) - (1{||C||<1}?, V5((x)) 



Xijflij 



< 



sj\{0}n{||^||<i} 



liePAl 

r f' d'g{x + sO ^^ ^^^^ 

j-^ \Jo dxkldXmn 



k,Lm,n 
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^s+\|o}n|i|fi|<i} m 



Js+\{o}n{m\>i} 

ll£l'^ 

's+\{o}n{i|cj|<i} 

+ iUx + Oi\\^ + m\\oo + \\9{x)x^J\U 

■J s+\{o}n{m\>i} 

xif^tj{dO + f^ij{dO) 

where Ci,C2 and C3 denote some constants and fJ-fj^Nj correspond to the 
Jordan decomposition fiij = fifj — . In the second last inequahty, we use 
the estimate Xij < ||x + ^||. The same as above can be shown for the measure 
m{dS,), whence (4.17) holds true. As by the first part of the proof, we have 
= A for all elements of A4, (4.17) implies 

lim \\Afn -A^f Woo = 0. 

n— >oo 

Since the infinitesimal generator of every Feller process is a closed operator, 
it foUows that S+ C V{A) and A = A^ on S+. □ 

4.3. Linear transformations and canonical representation. In this sub- 
section, we shall deal with linear transformations of afRne processes. The 
proposition below states how the parameters of an affine process on S'j' 
change under such linear maps, which allows us to establish a canonical 
representation of an affine process. 

Proposition 4.13. Suppose X is an affine process on S'^ with param- 
eters a,f5^^,c, 7,m,/u as specified in Definition 2.3 and b£ Sj". Further- 
more, let G : Sj" —7- Sj" , x 1— )• gxg~^ be an automorphism, where g G is in- 
vertible. Then, Y := gXg~^ is an affine process on , whose parameters, 
denoted by ~, are given as follows with respect to the truncation function 

x = gx{g~^i{9^)~^)9^ : 

b = gbg'^, 
c = c, 
m{d() = G^m{d$,), 
a = gag~^, 
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^^^^^ = (^^WFWt)^'''^''^*^^^^^'''' 

B^{u) = {g'r^B^{g'ug)g~\ 
where G^m (G^fi) is the pushforward of the measure m (n, resp.). 

Proof. Let us consider the process 

for which we have 

E[exp(-(n,y,^))] 

= E[exp{-{g''ug,Xf"y^'^"^"))] 

= exp{-(j){t,g~^ug) - {i;{t,g'^ ug),g~^y{g'^)~^)) 

= exp{-(l)(t,g'^ug) - {{g~^ y^ij{t,g~^ ug)g~^ ,y)). 

Define now (j) and ip by 

4>{t,u) = (l){t,g~^ug) and i){t,u) = {g~^)~'^'4){t,g~^ ug)g~'^ , 

to see that Y is an affine process on 5^. Using (2.14) and (2.16), we conse- 
quently obtain 

d4>{t,u) ^ d(t){t,g~^ug) 
dt ~ dt 

= F{i;{t,g'^ug)) 

= {b, ^{t, g~^ug)) +c- f (e-(V'(*.5^«9).5> _ l)m{di) 
Js+\{o} 

= {ghg'^Mt, u)) + c- [ (g- - l)m{dO 

Js+\{o} 

= {b,ip{t,u))+c- [ (e-^'^^*'")'^^ -l)G,m(dC). 
Js+\{0} 

Due to the uniqueness of the Levy-Khintchine decomposition, this imphes 
that b transforms to 6 = gbg~^ , c remains constant and m becomes fh[d^) = 
G^m{d£^). For -i/; we proceed similarly, that is, we have 
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= {g^)-^-2i;{t,g^ug)a^{t,g^ug) + {^{t,g'^ ug)) + 



.MoV MF^i J^^^'^T ' 

from which it can be seen that a transforms to 5 = gag~^ , 7 becomes 7 = 
{g~^ )~^jg~^ , and fi changes to 

for every E G S{S^ \ {0}). Moreover, since x = gx{9^^ ^{9^ )^^)9~^ 

(4.18) B^iu) = ig^r'B^ig'^ug)g~\ □ 

By means of Proposition 4.13, we can derive a canonical representation 
for affine processes. 

Proposition 4.14. Let X be an affine process on S'^ with parameters 
a,/3*-', c, 7, m,/i as specified in Definition 2.3 and b£ S'j'. Then there exists 
an automorphism G:S^ — )• S^,x i— )• gxg~^ such that the parameters of the 
affine process Y = gXg^ , denoted hy~, are as in Proposition 4-13 with 

b = e = diag{dn,---, Odd), a = lf, 



where we define 



Ir 




Proof. By Proposition 4.13, the parameters of y = gXg~^ transform as 

a = gag~^, b = gbg~^. 

Since a and b G S'j', they are jointly diagonalizable through an automor- 
phism on 5^" . More precisely, there exists an invertible matrix g G Md such 
that 

gag~^ = if with r = rk(Q;) 

and 

gbg'^ = diag(6'ii, . . . , Odd) =■ 0, 



where rk denotes the rank of a matrix. For the proof of this fact, we refer 
to [22], Theorem 8.7.1. □ 
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4.4. Condition on the constant drift. This subsection is devoted to show 
that condition (2.4) holds true for any affine process X on 5^". Since the 
automorphism G : Sj" — )• S'j' in Proposition 4.14 is order preserving, it suffices 
to consider affine processes of the canonical form as specified in Proposition 
4.14. The following result is a consequence of the Levy-Khintchine formula 
on R+. 



Lemma 4.15. Let Y be an affine process of canonical form as specified 
in Proposition 4-14 with parameters denoted hy~. Then, for any y G dS^ , 
we have 

(4.19) Vdet{y)eNs+{y), [ (x(0, V det(y))M(y, dO < oo 

Js+\{0} 

and 

(^,Vdet(y)) + (5(y),Vdet(y))- / (x(0, V det(y))M(y, 

Js+\{0} 

(4.20) ^ 

+ 2 X] ^ij ^ki det(y) > 0. 

i,j,k,l 

Proof. Let y £ dSj and let / G C^{S^) be a function with / > and 
f{x) = det(3;) for all x in a neighborhood of y. Then, for any v G IR+, the 
function x i— )• e~^^^^^ — 1 lies in C^(S'j') and thus in ^{A), where A denotes 
the infinitesimal generator of Y. Note that f{y) = 0. Hence, the limit 

^(e-^/(y) _ 1) = i / (e-''^^^) - l)pt{y,dO 
t^o+ t 



hm i / {e-'^'-l)p{{y,dz), 

^0+ t Jm 



t^o+ t 

exists for any v G M_|_, where Pt{y,d(,) denotes the transition function of Y, 
and p{{y,dz) = f^pt{y,dz) is the pushforward of pt{y, •) under /, which is a 
probability measure supported on M+. 

Using the same arguments as in Proposition 4.9 [i.e., applying Lemma 4.5 
as done below (4.4)], and noting that f{y) = 0, we conclude that 

V ^ A{e-''f'^y^ - 1) 



(4.21) = ^ E ^ijki{y){v^di,f{y)dkif{y)-vdi,dkifiy)) 

i,j,k,l 

-v{e + B{y),Vf{y)) 
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+ / (e--^(^+«)-l + 7;(x(e),V/(y)))M(2/,dC) 
Js+\{o} 

is the logarithm of the Laplace transform of an infinitely divisible distribu- 
tion on M4.. Note that 



(Vdet(y),2;) = ^det(y + ix) 



> 0, xG 5+, 
t=Q( = 0, x = y. 



Hence, Vdet(y) G Ng+{y) and the admissibility condition (2.9) implies (4.19). 

By the Levy-Khintchine formula on M-|_ (see [49], Theorem 3.21), the linear 
coefficient in v in (4.21) has to be nonpositive. But this is now just (4.20), 
whence the lemma is proved. □ 

It now remains to show that (2.4) follows from (4.20). For this purpose, 
it suffices to evaluate (4.20) at diagonal elements y S dS'^ . Thus, we state 
the following lemma. 

Lemma 4.16. Let y e Si be diagonal, and let f G C^{Sf ). Then we 



have 



1 

- (yik ilr)jl + yu {lt)jk + Vjk {Ir)il + Vjl 



i,j,k,l=l 



1 



dx% 



ik 



+ 



dxij dxki 
d'fix) 



x=y 



x=y dXijXji 



Proof. Obvious. □ 

Next, we calculate the partial derivatives of the determinant 



x=y 



Lemma 4.17. Let y £ be diagonal, y = d\a.g{yi\,y22, ■ ■ ■ ,ydd)- Then 
we have 

Wvkk, 

k^i 

0, 



dAei{x) 



dx 



x=y 



else. 



and 



52 det(x) 



dxijXji 



92 det(a 



x=y k=l,k^i,k^j 



dx% 



for I < i < j < d, 
forl<i<j<d, 



x=y 
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where the empty product is defined to be 1. 

Proof. In dimension d = 2, the assertion is easily checked, as det{y) = 
2/112/22 - 2/122/21- Therefore, we have 

(9iidet(y) =2/22, <922 det(y) = yn, <9i2 det(y) = ^21 det(y) = 

as well as 

dj, det(y) = a|2 det(y) = df^ det(y) = det(y) = 0, 

di2 821 det(y) = ^21 9i2 det{y) = -1. 

For dimension strictly larger than 2, we employ a combinatorial argument. 
Recall Leibniz's definition of the determinant. 



(4.22) 



det(x) = ^ sgn(cr) 



0-6S 



k=l 



where a is an element of the permutation group S on the set {1, 2, . . . ,d} 
and sgn denotes the signum function on S, that is, sgn = 1 if u is an even 
permutation and sgn= —1 if it is odd. Differentiation of (4.22) with respect 
to Xij yields 



ddet{x) 



dxi 



x=y 



J^sgn((T)l|,(,)=j| 



Wvkk, 



x=y 



k^i 

0, 



Thus, for the second derivative we have 
d^det{x)_ 



x=y 



^sgn(a)l|^(,)=j}l|^(j)=i| Yl Xfc^(fc)j 



k^i^j 



x=y 



else. 



kj^ij^j 



where the last equality holds since y is diagonal. For dfj det(x), the statement 
is obvious. □ 

We are prepared to prove the admissibility condition on the constant drift. 



Proposition 4.18. Let X be an affine process on then (2.4) holds, 
that is, 

bt{d-l)a. 

Proof. Since the automorphism GiS'^ — t- S'^ in Proposition 4.14 is 
order preserving, it suffices to show that (4.20) in Lemma 4.15 implies 

(4.23) ey{d-i)i^. 
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We show that 6mm ^d — 1, if r>m. To this end, take again some diagonal 
y £ dSj of form y = diag(yii > 0, . . . , ymm = 0, . . . , y^rf > 0). By Lemmas 4.16 
and 4.17, we obtain 



^ 9ii da det(y) + Y,iB{y))ij ^^j det(y) 



1=1 



^Y,{x{0)ijdijdet{y))M{y,dO 



+ ^ ((yiil{i<r-} + yij'^{i<T}){Q% det(y) + 5^ dji det(y))) 
= [Qh n Vkk\ + {fLmVn n Vkk) 

sr^ f {x{0)m7nyilY\k^mykk ^ , 

— BFx^ — 

" ^ ( n yfcfcl{i<r} + n ykkhi<r} ) 
i^j k^i ^ 

= 9mm \\ykk+\\ ykk (V (Ptmyu " Vu [ ^^Stt^w; (f^C) ) ) 

k^m k^m Vm^ ^5+\{0} 11^ II A 1 7 7 

- (^^ - 1) n ykk'^{m<r} > 0- 

As YlkjtmVkk > and by (2.11) also 

i^^myu-yu i WTTMdo)>^ 
\ Js+\{a} ll?lrAi 7 

for I / m, letting ya ^ 0,1 ^ m yields 8mm > d — 1 ior r > m. Relabeling of 
indices then proves (4.23). □ 

5. Sufficient conditions for the existence and uniqueness of affine pro- 
cesses. In this section, we prove that for a given admissible parameter set 
a, b, f]"^^ , c, 7, m, ji satisfying the conditions of Definition 2.3, there ex- 
ists a unique affine process on Sj", whose infinitesimal generator A is of 
form (2.12). Our approach to derive this result is to consider the martingale 
problem for the operator A. In order to prove uniqueness for this martingale 
problem, we shall need the following existence and uniqueness result for the 
generalized Riccati differential equations (2.14) and (2.15). 



38 CUCHIERO, FILIPOVIC, MAYERHOFER AND TEICHMANN 

5.1. Generalized Riccati differential equations. We first derive some prop- 
erties of tlie function R given in (2.17). 

Lemma 5.1. R is analytic on S'^^ and quasi-monotone increasing on 
"^d ■ 

Proof. That R is analytic on Sj"^ follows by dominated convergence 
(see, e.g., [16], Lemma A. 2). 
Now let 5 > 0, and define 

R'{u) = -2uau + B^{u)+^- [ (£:^!l_I±MM V(dO 



-2uau + 7 + ^-B^(ii) - j 



Now, the map u i— )■ —2uau + 7 is quasi-monotone increasing, as it is shown 
in Step 3 of the proof of Proposition 4.9. Furthermore, it follows from the 
admissibility condition (2.11) that 



J\m\>&\ llClr Al 



is a quasi- monotone increasing linear map on S'j" . Finally, the quasi-monotonicity 
of 

is a consequence of the monotonicity of the exponential and that supp(/.f) C 
"^d ■ 

By dominated convergence, we have \\m.^^QR^[u) = R{u) pointwise for 
each u S St . Hence, the quasi-monotonicity carries over to R. Lideed, choose 
x,u,v G such that u ^ v and {v — u,x) = 0. Then we have for all 6, 
{R\v)-R\u), x) > 0. Thus, 

{R\v) - R\u),x) {R{v) - R{u),x) > 0, 

as 5 — 7- 0, which proves that R is quasi- monotone increasing. □ 



Lemma 5.2. There exists a constant K such that 

K 
~2 



(5.1) (^,i?(^))<_(||^||2 + i), ^^sj. 
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Proof. We may assume, without loss of generality, that the truncation 
function in Definition 2.3 takes the form xiO — ^{U\\<i}^ [otherwise adjust 
B{u) accordingly]. Then, for all u G Sj" we have 

Riu) = -2uau + B^iu)+j- [ ('£:!!l_i±iil!^) 

■''5+\{0}n{|lc|l<i} V 114 II 



>o 



(5.2) 



•''s+\{o}n{i|5!l>i} 



^ -2uau + 5 ' (n) + 7 + /i(5+ n {||^|| > 1}) 

where we use that 

5;\{0}n{||c|l>i} ^s+\{o}n{||^||>i} 
Set now 

7:=7 + M5in{||e||>l})G5+. 

By (5.2), we obtain, for u G Sj", that 

{u,R{u))<{u,B^{u)) + {u,j), 

from which we derive the existence of a positive constant K such equation 
(5.1) holds. □ 

Here is our main existence and uniqueness result for the generalized Ric- 
cati differential equations (2.14) and (2.15). 

Proposition 5.3. For every u G S'^'^ , there exists a unique global M4. x 
S^~^ -valued solution of (2.14) and (2.15). Moreover, 4){t,u) andtlj{t,u) 

are analytic in {t,u) G x S*^"*". 

Proof. We only have to show that, for every u G 5'j"'', there exists a 
unique global 5j"''-valued solution ip of (2.15), as then (/) is uniquely deter- 
mined by integrating (2.14) and has the desired properties by admissibility 
of the parameter set. 

Let u G •S'j"''. Since R is analytic on S'^^ , standard ODE results (e.g., 
[14], Theorem 10.4.5) yield there exists a unique local 5j"''-valued solution 
'tpit,u) of (2.15) for t£ [0,i+(u)), where 

t+(n) =lim inf {t > | \\'ip{t,u)\\ > n or Tp{t,u) G dSj} < 00. 
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It thus remains to show that t^{u) = oo. That i{j{t,u), and hence (f){t,u), is 
analytic in (t,n) G M+ x then follows from [14], Theorem 10.8.2. 

Since R may not be Lipschitz continuous at dS^ (see Remark 5.4 below), 
we first have to regularize it. We thus define 

R{u) = -2uau + B''{u)+^- [ /'^J!l_l±iii^ ^(^0- 

■/s+\{o}n{||5||<i}V ll^lr J 

It then follows as in Lemmas 5.1 and 5.2 that R js quasi-monotone increasing 
on Sj" and that (5.1) holds for some constant K. Moreover, R is analytic on 

Sd- Hence, for all u G 5"^, there exists a unique local S'rf-valued solution ijj of 

^^^ = R(4,{t,u)), ^(0,n)=^, 

for all t G [0,t^{u)) with maximal lifetime 

t+(n)=lim inf {t>0 \ \\ij{t,u)\\ >n}<oo. 

n— >oo 

From (5.1), we infer that for all u G and t < t^{u), 

at||V^(t,n)||2 = 2(^(t,n),9t^(t,n)) <^(||^(t,^z)f + 1). 

Gronwall's inequality (e.g., [14], (10.5.1.3)) implies 

(5.3) \m,u)f<e''\\\uf + l), t<I+iu). 

Hence, t^{u) = oo for u G S'^ . As R is quasi- monotone increasing on Sj", 
Volkmann's comparison Theorem 4.8 now yields 

Q<ij{t,u)<^{t,v), t>0, for allO^u^v. 

Therefore and since ij){t,u) is also analytic in u. Lemma 3.3 implies that 
i^{t,u) G Sj^ for all (t,n) G M+ x 5++. 

We now carry this over to ip{t,u) and assume without loss of generality, 
as in the proof of Lemma 5.2, that the truncation function in Definition 2.3 
takes the form xiO = 1{||5||<1}'^- Then 

R{u) - R{u) = - [ (e-<"'«> - l)fi{dO to, u G 5+. 

Js+\{o}n{m\>i} 

Hence, for u G S'j"*' and t < i+(tt), we have 

- mt.u)) = ^ - Rm.u)) :< ^ - Rm.u)). 

Theorem 4.8 thus implies 

iP{t,u)yi^{t,u)eS++, tG[0,i+(n)). 
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Hence, t+{u) = liminf„_^.oo{i > | [['(/'(ijit)!! > n}. Using (5.1) again, we now 
can show as for ip that 

||^(t,n)|p<e^*(||u||2 + l), t<Uiu). 

Hence t^{u) = oo, as desired. □ 

Remark 5.4. Lemma 5.1 states that the admissibihty of the parame- 
ters a,/?*-', imphes quasi- monotonicity of R on S'j'.^ Moreover, quasi- 
monotonicity just means that R is "inward pointing" close to the boundary 
S'j". Indeed, let u,x ^ 5'j' with {u, x) = 0. Then {R{u),x) > (7, x) > 0. Hence, 
if R were Lipschitz continuous on , a deterministic variant of Theorem 
A. 5 would imply the invariance of Sj" with respect to (2.15) right away. How- 
ever, the map R might fail to be Lipschitz at dS'^ (see the one-dimensional 
counterexample [16], Example 9.3), even though it is analytic on the inte- 
rior iSj"*". Here, quasi-monotonicity plays the decisive role. It leads to the 
phenomenon that 'ip{t,u) stays away from the boundary dS'^ for u S 5*^"^, 
which is of crucial importance in our analysis. 

5.2. The martingale problem for A. We are now prepared to study the 
martingale problem for the operator A given by (2.12). For the notion of 
martingale problems, we refer to [17], Chapter 4. We shall proceed in four 
steps. First, we approximate A by regular operators A'^'^'"' on the space 5+ 
of rapidly decreasing C°°-functions on S'j", defined in (B.l). Second, using 
Theorem A. 5 below, we show that there exists an S'j'-valued cadlag solution 
of the martingale problem for A^'^'"'. Third, a subsequence of these solutions 
is shown to converge to an Sj" U { A}-valued cadlag solution of the martingale 
problem for A. Finally, we show that this solution is unique, Markov and 
affine, as desired. 

Note that we cannot employ Stroock's [51] seminal existence and unique- 
ness results for martingale problems, since those are solved on and require 
uniform elliptic diffusion parts. Neither of these is satisfied in our case. 

Now let (a, 6, /S*-' , c = 0, 7 = 0, m, fi) be some admissible parameter set. Fix 
some e,6 > and n G N. In order to bound the coefficients and cut off the 
small jumps, we let 

<^„GCr(5d), 

(5.4) 

^We conjecture that the converse also holds: R is quasi-monotone on Sj" and supp(^) C 
Sj" if and only if the parameters a, 7, are admissible. 



0<iPn<l, 




\x\\ < n, 
|x|| > n + 1. 



42 CUCHIERO, FILIPOVIC, MAYERHOFER AND TEICHMANN 

We then define the bounded and smooth parameters 

m'^(dC) = "T'('^Ol{||5||>5}, 

Concerning the diffusion function Aijki{x) given by (2.13), we first find 
an appropriate factorization which will allow us to write the continuous 
martingale part of X as a stochastic integral. Thereto observe that any Sj"- 
valued solution, presumed that it exists, of the following symmetric matrix- 
valued diffusion SDE: 

(5.5) dZt = y%dWt^ + ^'^ dW^ 

where is a standard d x d-matrix Brownian motion and S G M^, with 
= a, has quadratic variation d{Zij, Zki)t = ^ijki{Zt)- Define now a^\x) G 
Sd by 

(5.6) a^\x) = V^M'^'S + yJm^^^, 
where M^^l=6ik6ji. Then (5.5) can be written as 

d 

dZt= <y'''{Zt)dWt,ki 

k,l=l 

and = Eln=i^r(^)^r(^)- 

Since a^\x) involves the matrix square root, which is neither Lipschitz 
continuous nor bounded nor globally defined, we need to introduce some 
approximating regularization in order to meet the assumptions of Theorem 
A. 5 below. Thereto fix some truncation function 



and define 

(5.7) Se,n{x) 



r]s{ipn{x)x){y^ipn{x)x + eld - Veld), if x e - eld 
0, otherwise. 



Note that Ss^n satisfies: 

• Ss,n G C^{Sd, Sd), 



Se,n{x) = {^y(pn{x)x + eld - VeT2) OU Sj, 
limj^O+ S£,n{x) = y^ipnix)x. 
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With this, we can now define the regularization of cj'^' by 

(5.8) af^,X^) = Se,n{x)M^^i: + STM''=Se,n(x), 

which then satisfies the smoothness condition of Theorem A. 5. Finally, we 
set 

d 

m,n 

(5.9) = {sl^{x))ikaji + {sln{x))iiajk 

+ (Se,n(a^))jfcaii + {sl^{x))jiaik, 
and define the corresponding regularized operator on Co{Sd) 



+ Y^{h,^ + B:^{x)) 



df{x) 



dxi 

(5.10) 



/ {f{x + i)-f{x))m'{di) 

Js+\m 



+ 

+ I (/(x + 0-/(x)-(x(e),V/(a;)))M^'"(x,dO. 

We now show that A^'^'"^ approximates A. We let S = S{Sd) and 5+ 
denote the locally convex spaces of rapidly decreasing C°°-functions on Sd 
and 5j" defined in (B.l) below, respectively. 

Lemma 5.5. S C D(^^''''") and, for every f £S+, 
(5.11) lim||^^'^'"/-^/||oo = 0. 

e,o,n 

Proof. Since (fn as defined in (5.4) converges uniformly on compact 
sets to 1, this is clear for the differential operator part. Concerning the 
integral part, we have 

/ (/(x + - fix) - (l|||5,|<i}^, V/(x)))(M'5'"(x,(ie) - M(x,dO) 
Js+\{0} 



< 



/(x + 0-/(x)-(l{||5||<i}e,V/(:E)) 



Ir^nmV lief A 1 
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+ 



E 



s7\{o} 



/(x + e)-/(x)-(l|||g||<i|g,V/(x)) 
II^PAl 



By dominated convergence the second term goes uniformly in x to 0, thus 
we only have to consider the first one. By splitting the first integral into 
■All€ll<i} ~'~ -All^l^i}' ^ote that || /{j|^||<i} II converges uniformly in x to 0. 
Hence, it remains to analyze 



E / {f{x + ^)-f{x))xi,{Mx)-l)fii,{dO 



which can be estimated by 



E 



{||CII>i} 



\\f{x + ^)x,,{^n{x) - l)\MidO+f^7M)) 



+ 



{||CII>i} 



\\f{x)x,,{y,n{x) - imf^tM)+f^^M))), 



where fifj , correspond to the Jordan decomposition of fiij = l^-tj ~ l^ij ■ 
/ lies in 5+, the second term converges uniformly to 0. For the first one, 
observe that for every n 

ii/(x + - 1)11 < \\f{x + o^ijW < \\f{x + m Ik + eii, 

such that we can apply dominated convergence. Again, since / lies in 5+, 
the first integral converges uniformly in x to as well. Hence (5.11) holds 
true, and S C 'D{A'''^'"') follows similarly. □ 

We now establish existence for the martingale problem for A'''^'"'. 

Lemma 5.6. For every x G 5'j' there exists an S'^ -valued cadlag solution 
X to the martingale problem for A^'^'"^ with Xq = x. That is, 

f{Xt)- fA''''^f{Xs)ds 
Jo 

is a martingale, for all f £ S. 
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Proof. Consider the following SDE of type (A.l): 



(5.12) 



Xf'"" = x+ f(h + B'\Xf''') - [ xiOM^'^'iXf''', ds 
Jo \ Js+\{o} J 



TT Jo 



k,l 

where W is a d x d-matrix of standard Brownian motions and J a finite 
activity jump process with compensator m^{d£,)+M^''^{X^'^'"', d^). Note that 
the quadratic variation of the continuous martingale part of X^'^'"' is given 
by A^^^i{x) as defined in (5.9). It thus follows by inspection that any cadlag 

solution X^'^'"' of (5.12) solves the martingale problem for A"^'^'". 

Hence, it remains to show that there exists an S'j'-valued cadlag solu- 
tion of (5.12). Let us recall the normal cone (2.23) to Sj". As 6 + B^{x) — 
fs+\{o}X{0 X M'5'"(x,(iO, o-^f„(x) and m^{dO + M'5'"(x,dO are designed 

to satisfy the assumptions of Theorem A. 5 and since x + supp(m''(-) + 
M^'"'(x, ■)) Q for all x G S^, we only have to show that for all x G dS^ 
and u £ Ng+ (x) 

(5.13) {af^,,{^),u)=0, 

/b + B^{x)- [ x{OM''''{x,dO 
, , \ Jsl\{0} 

(5.14) 

-^E^^'!n(^Kn(x),^)>0. 

k,l=l I 

Due to the definition of ^"^'^(x), respectively, the definition of Se^n{x) given 
in (5.7), condition (5.13) is satisfied. Concerning (5.14), we have by (2.11) 

B^{x)- I xim''''i^,dO,u)>0. 
Js+\{o} I 

Moreover, it is shown in Lemma 5.7 below that 

{h-\Y.Daf^^{x)al^{x\^>^. 

\ k,l=l I 

The lemma now follows from Theorem A. 5. □ 

Lemma 5.7. Let x = OKO^ G S'^ where A = diag(Ai, . . . , \d) contains 
the eigenvalues in decreasing order and let o"^'^ be defined by (5.8). Then, 
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for all X G S][ , 



^ly^ 99n(x)(v^Ajy?n(x) +g- Ve) 

2 ^ ^ Ai99„(x) +£ + ^J XjiPn{x) + E 



2 ^ 2Y^Ai99„(x) + e 
w/iere = ((STS)0)„,0„i + ((STs)0)„iO™ and Z^^ = O^OfciE;, 

OniOlsiYlm ■ 

Furthermore, if 

(5.16) 5^(d-l)S^S, 

(5.17) ^ E ^<n(^Xn,(^),^\ > 
\ k,l=l I 

for all X G dS'^ and for all u € Ng+ {x) . 
Proof. Let us denote 

1 



k,l=l 

and notice that 



ki 



We now use the following formula from [29], Theorem 6.6.30: 



f{V{t)) = 0{t)(yAf{X,{t),\,{t))AnO{tyv'{t)0{t)^^^^^ 



dt 



where V(t) = 0{t) diag(Ai(t), . . . , \d{t))0{t)~^ is a family of symmetric ma- 
trices and Af{u,v) = for u^v and A/(n,n) = f'{u). This holds 

true if V{-) is continuously differentiable for t G {a,b) and /(•) is continu- 
ously differentiable on an open real interval which contains all eigenvalues 
of V{t) for ah te{a,b). 
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We now apply this formula to our case, where f(t) = \/t and 

V{t) = ^n{x + taf^n{x)){x + taf^nix)) + sh- 

Since we take the derivative at t = 0, we only have to consider 

V{Q) = 0{^n{x)K + eh)0^, 
where O is the orthogonal matrix diagonalizing x and 

Note that we do not have an explicit contribution of % which is part of the 
definition of since r?e(5j') = 1 and Vr/£(5'j") = 0. Some lines of calcula- 
tions then yield (5.15). 

Let us now verify (5.17). Take an arbitrary x = OKO^ G dS'^ and assume 
first that it has rank d — 1, that is, = and all other eigenvalues are 
strictly positive. By Lemma 4.1 and (2.23), the elements of Ng+{x) can 

then be written as u = OKO^ , where K = diag(0, . . . , 0, /c) with > 0. Thus, 
(5.17) now reads 

{h - C''~^{x),OKO'^) = kp'^hO - O^C"'"(x)0]dd. 

As [0^^7*0]^^ = 26id{0'^i:^^0)id and Z^^'O = 26idOti{^0)id, we have 



Since V, , . ^fni^K^^p+i^ < _ l, we obtain by condition (5.16) 

[O'^bO - 0'^C''"{x)0]dd > [O'^ib -{d- 1)STS)0]^^ > 0, 

which proves (5.17) for x G dS^ with rk = d — 1. In the general case, we can 
proceed similarly. For x G dS'^ with rk = r < d — 1, the elements of 



X 



are given by n = OKO^ , where 

K- 



d 




k 



with k G S'^_^.. This follows again from Lemma 4.1 and (2.23). Now, (5.17) 
can be written as 

(6-C"'"(2;),OKO^) 



> (0^(6-rS^S)0,K) >0, 
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which proves the assertion. □ 

Combining Lemmas 5.5 and 5.6, we obtain the announced existence result 
for the martingale problem for A. 

Lemma 5.8. For every x G S'j', there exists an S'j' U {A} -valued cadlag 
solution X to the martingale problem for A with Xq = x. That is, 

f{Xt)- f Af{Xs)ds 
Jo 

is a martingale, for all f G 5+ . 

:,<5j 



Proof. By Lemma 5.6, there exists a solution to the martingale 

problem for A'''^'"' with sample paths in ©(Sj") (the space of 5j'-valued 
cadlag paths), and hence also in ©(tSj" U {A}). We now claim that (X^-'^-'^) 
is relatively compact considered as a sequence of processes with sample 
paths in ©(Sj" U {A}).^'^ For the proof of this assertion, we shall make use of 
Theorems 9.1 and 9.4 in Chapter 3 of [17]. In order to meet the assumption 
of [17], Chapter 3, Theorem 9.4, we take C^{S'^) as subalgebra of Cfc(S'j'). 
Then, for every T > and / G (7^(5^"), we have 



sup 'Ex 

E,5,n 



essup|^^'^'"/(Xt' 
■te[o,r] 



< oo, 



since there exists a constant C such that ||^'^''^'"/||oo < Cp3,+{f) < co for all 
n,e,5, where pk,+ are the semi-norms as defined in (B.2) (see also the proof 
of Proposition 4.12). Thus, the requirements of [17], Chapter 3, Theorem 
9.4, are satisfied. Note that Y in the notation of [17], Chapter 3, Theorem 
9.4, corresponds in our case to f{X) such that [17], Chapter 3, Condition 
(9.17), is automatically fulfilled. It then follows by the conclusion of [17], 
Chapter 3, Theorem 9.4, that (/(X|'^'")) is relatively compact [as family of 
processes with sample paths in D(M)] for each / G C^(S'j'). Furthermore, 
since we consider Sj" U {A}, the compact containment condition is always 
satisfied, that is, for every r] > and T > 0, there exists a compact set 
r^,T C (5+ U {A}) for which 

inf Fx[X'/''' G Te,T for t G [0, T]] > 1 - ry 

e,S,n ' 



^"This means that the family of probabihty distributions associated to (X^'*'") is rel- 
atively compact, that is, the closure of (P^'^'") in 7'(D(S'j" U {A})) is compact. Here, 
P(D(5j" U {A})) denotes the family of probability distributions on D(Sj U {A}) and 
the distribution of X""'^'". 
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holds true. By [17], Chapter 3, Theorem 9.1, and the fact that {1, C^(5j")} is 
dense in C(5j'U{A}), we therefore obtain that (X"^'^'*^) is relatively compact 
in ^{S^ U {A}). Thus, there exists a subsequence (p^fe'^fe'"*) of the proba- 
bihty distributions associated to (X*"'^'") which converges in the Prohorov 
metric to some hmit probabihty distribution. By [17], Chapter 3, Theorem 
3.1, this implies weak convergence of (P'^fc'^fc'"fe) and hence the subsequence 
(X^'"'^'""'*) converges in distribution to some hmit process X in B(S'j'U{A}). 

Combining this with Lemma 5.5 and [17], Chapter 4, Lemma 5.1, we 
conclude that X is a sohition to the martingale problem for A. Hence, the 
lemma is proved. □ 

We can now prove the existence and uniqueness of an affine process for 
any admissible parameter set. 

Proposition 5.9. Let {a,b, 13^^ ,c,^,m, fi) be an admissible parameter 
set. Then there exists a unique affine process on S'^ with infinitesimal gen- 
erator (2.12), and (2.1) holds for all {t,u) G M+ x S'^ , where (j){t,u) and 
ip{t,u) are given by (2.14) and (2.15). 

Proof. Suppose first that c = and 7 = 0. Let x G S'j'. Then Lemma 
5.8 implies the existence of an 5^ U {A}-valued cadlag solution X of the 
martingale problem for A with Xq = x. We now show that X is unique in 
distribution. 

Thereto, note that by [17], Chapter 4, Theorem 7.1, 

(5.18) f{t,Xt)- [\Af{s,Xs) + dsf{s,Xs))ds 

Jo 

is a martingale for all rapidly decreasing functions / G 5(M4. x /Sj"), similarly 
defined as 5+ in (B.l). Now let 4> and ip be the unique solutions of the gener- 
alized Riccati differential equations (2.14) and (2.15), given by Proposition 
5.3. Fix t > 0, ti G Sj~^, and some / G 5(R+ x Sj) such that 

/(s, x) = e-'^(*-^.«)-('/'(*-^.«).^) , 0<s<t,xeS+. 

Then 

Af{s,x) + dsf{s,x) = 0, 0<s<t,x€ Sj. 
In view of (5.18), the Laplace transform of Xt at u is thus given by 

(5.19) E^[e-<"'^'>] =E^[/(t,Xt)] = /(0,x) -0 = e-'^(*'")-<^'(*'")'^>. 

Since u G 5'j"'' was arbitrary, we conclude that the distribution of Xt is 
uniquely determined for all t > 0. From [17], Chapter 4, Theorem 4.2, we 
infer that X is a Markov process with generator A on and thus unique in 
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law as solution of the martingale problem for A. Moreover, by (5.19), X is 
stochastically continuous and affine. Thus, the proposition is proved under 
the premise that c = and 7 = 0. 

For general parameters c and 7, we employ a Feynman-Kac argument. 
Denote by B and (Qt) the affine generator and corresponding Feller semi- 
group associated with (a, b, , c = 0, 7 = 0, m, ^u) from the first part of the 
proof, respectively. Since x 1— )■ c+ (7, x) is nonnegative on S'^ , it follows along 
the lines of [16], Proposition 11.1, that 

defines a Feller semigroup (Pt) on Co(5'j') with infinitesimal generator Af{x) = 
Bf{x) — (c+ (7, for / E 5+, which is the desired solution. □ 

5.3. An alternative existence proof for jump processes. For affine pro- 
cesses without diffusion component (i.e., the admissible parameter a van- 
ishes), the existence question can be handled entirely as in the case of affine 
processes on WJ^ x [16], Section 7. In this section, we elaborate an al- 
ternative existence proof in this specific case, by following the lines of [16]. 
Note that the OU-type processes driven by matrix Levy subordinators [3] 
are contained in the class of pure jump processes of this section. 

We call a function / : S'j' — > M of Levy-Khintchine form on , if 

f{u) = {bo,u)- [ (e-<"'«> - l)mo(dO, 
Js+\{0} 

where 60 £ S'^ and mg is a Borel measure supported on S'j" such that 

/ (IICII Al)mo(dO<oo. 
Js+\{0} 

Once again, we recall that a distribution on Sj" is infinitely divisible if 
and only if its Laplace transform takes the form e~-^^'^\ where / is of the 
above form (see also Step 1 in the proof of Proposition 4.9). 

Similarly to [16], we introduce the sets 

C := {/ + c I / : 5'j' — t- R is of Levy-Khintchine form on tSj", c G 

Cs '■= {ip \ {Tp{u),x) € C for all x G "Sj"}. 

The following technical statement can be obtained easily by mimicking the 
proofs of the corresponding statements in [16], Proposition 7.2 and Lemma 
7.5: 

Lemma 5.10. We have: 
(i) C, Cs are convex cones in C(S'j'). 
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(ii) (/> G C, t/j ^Cs imply cl){ip) G C. 

(iii) tp,ipi £ Cs imply ^pl{^p) GC5. 

(iv) If (pj^ £ C converges to a continuous function cj) on S'^ , then (j) ^C. 
A similar statement holds for sequences in Cs- 

(v) Let (a = 0, c, 7, m, /i) be an admissible parameter set. Then 
^ R locally uniformly as (5 — )• 0, where R^ corresponds to the admis- 
sible parameter set (a = 0, c, 7, m,/il|||^||>5}). (Note that there is one 
fixed truncation function.) 

Proposition 5.11. Let {a = 0,b, /3^^ ,c,'y,m, fi) be an admissible param- 
eter set. Then for all t > , the solutions {(f){t,-),'ip{t,-)) of (2.14) and (2.15) 
lie in {C,Cs). 



Proof. Suppose first that 



11 



(5.20) / < oo 



/5+UO} llell Al 

for all i <j. Then equation (2.15) is equivalent to the integral equation 



ft _ 

(5.21) ^P{t,u) = e^^\u)+ / e-^' (*-^)i?(^(s, u)) ds, 

Jo 

where R{u) = R{u) + {u) and G C{Sd) is given by 
B\u):=B\u)- [ M|M;,(de). 

Js+\{o} IICIrAi 

Here, e^^^{u) is the notation for the semi-group induced by dtx{t,u) = 
B'^{x{t, u)), x{0 ,u) = u. Hence, the variation of constants formula yields 

(5.21) . 

Due to admissibility condition (2.11), we have that B~^ is a linear drift 
which is "inward pointing" at the boundary of Sj", which is equivalent to 

c^^* being a positive semi-group, that is, e^^* maps 5j" into S^. Therefore, 

e^^^ G Cs and since R{u) is given by 

Js+\{o} IkII Al 

with /i satisfying (5.20), we also have 

(5.22) R^Cs. 



According to our conjecture in Section 2.1.4, this would already cover all possible 
jump measures if d > 2. 
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Using Picard's iteration and Lemma 5.10, it follows that the sequence 
t/^^'^^ defined as 

■.= u, 

Jo 

lies in C5, for each t > 0, hence so does its limit Tp{t,-). Since F gC, we have 
agam by Lemma 5.10 (j){t, •) = £ ^(^'(s, •)) ds G C. 

By an application of Lemma 5.10(v), the general case is then reduced to 
the former, since clearly satisfies (5.20). □ 

We are prepared to provide an alternative proof of Proposition 5.9 un- 
der the additional assumption a = 0: by Proposition 5.11, {(l){t, ■),ip{t, ■)) 
lie in {C,Cs)- Hence for all t > 0, x G Sj", there exists an infinitely divisible 
sub-stochastic kernel ^^(x, d^) with Laplace-transform e"''^^*''")"^'^^*'")'^^ . The 
Chapman-Kolmogorov equations hold in view of properties (3.1) and (3.2). 
Whence, Proposition 5.9 follows. 

Remark 5.12. We note that the proof of statement (v) in Lemma 5.10 is 
much easier than the one of [16], Lemma 7.5, because a = 0. However, for a 7^ 
and d>2,R cannot be locally uniformly approximated by functions R^ of 
a pure jump type such as in Lemma 5.10. Indeed, otherwise one could infer 
as above the existence of an affine process which is infinitely decomposable 
and has nonvanishing diffusion component. This is in contradiction with 
Proposition 2.9 and in the case of pure diffusions it contradicts Example 
2.8. 

6. Proof of the main results. 

6.1. Proof of Theorem 2.4- The first part is a summary of Propositions 
3.4, 4.9, 4.12 and 4.18. The second part follows from Proposition 5.9. 

6.2. Proof of Theorem 2.6. Let X be a conservative affine process. It is 
shown in Proposition 4.12 that {e~^"''^ | u S S'j"''} C D{A). Hence, 

Jo 

is a (J-(,PxO -martingale with J^t defined in (2.18). From [31], Theorem II.2.42, 
combined with (4.16) and Remark 2.5, it then follows that X is a semi- 
martingale with characteristics (2.19)-(2.21). The canonical semimartingale 
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representation ([31], Theorem II. 2. 34) of X is thus given by 

Xt = x + Bt + Xf + f [ x{Oil^''{ds,dO-Hds,dO) 
Jo Js+\{0} 



+ f I {C-X{0)f^''{ds,d0, 
Jo Js+\^o} 



lo JS^\{0} 

where X'^ denotes the continuous martingale part, and the random mea- 
sure associated with the jumps of X. In order to establish representation 
(2.22), we find it convenient to consider the vectorization, vec(X'^) S , 
of X'^. The aim is now to find a d^-dimensional Brownian motion on a 
possibly enlarged probability space and a x d^-matrix-valued function a 
such that 

(6.1) vec(Xf)= f a{X,)dWs. 

Jo 

Thus, a has to fulfih 

d{Xfj,X^i)t = Xt^ikaji + Xt^iiajk + XtjkOiu + Xtjiaik 



(6.2) 



{a{Xt)a' {Xt))ijki. 



As suggested by (5.6), we define the entries of the d'^ x d^-matrix (j{x) in 
terms of a^\x) given in (5.6) by 

(6.3) Tijkiix) = o-fj(x) = Vxik^ij + ^JiVxjf,. 

Note that the (A;/)th column of cr{x) is just the vectorization of the ma- 
trix a'^^x). We thus obtain Aijki{x) = {a{x)a~^ {x))ijki- Hence, a{x) satisfies 
(6.2). Analogous to the proof of [45], Theorem 20.1, we can now build a 
d^-dimensional Brownian motion W on an enlargement of the probability 
space such that (6.1) holds true. As the (ij)th entry of X'^ is given by 

pt d ^ 
^t,ij = vec(Xf )ij = / V (JijkiiXs) dWsM 
■^0 k,l=l 

= [ {^/x'sdWs^ + ^'^ dWJ^/X, 
Jo 



where W is the d x d-matrix Brownian motion satisfying vec(iy) = W, we 
obtain the desired representation. 
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6.3. Proof of Theorem 2.9. We first prove some technical lemmas. 

Lemma 6.1. Let g-S^ — )■ M be an additive function, that is, g satisfies 
Cauchy 's functional equation 

(6.4) g{x + y)=g{x)+g{y), x,yeS^. 

Then g can he extended to an additive function f : Sd Moreover, if g is 
measurable on S'j' then f is measurable on Sd. In that case, f is a continuous 
linear functional, that is, f{x) = {c,x) for some cG Sd- 

Proof. The first part follows from Lemma 4.4. 

Concerning measurability, let E G i3(M), a Borel measurable set. Then we 
have by the additivity of /, 

oo oo 

f-\E)=\J Bn=[J{x + nId\xeSdJ{x)eE,\\x\\<n}-nId 

n=l n=l 

oo 

= (j{yeSd\ f{y) eE + f{nld), \\y - nhW < n} - nh 

n=l 

oo 

= [J{y G Sj" I g{y) e E + g{nId), \\y - nld\\ <n}- nid, 

n=l 

which is again a measurable set, in view of the measurability of g on S'j'. 

For X £ Sd we write x = {xi)i, where 1 < i < i^i^^tii. w/g introduce the 
additive functions /j :M — >• M via fi{xi) = /(O, . . . , 0, Xj, 0, . . . , 0). By the just 
proved measurability of /, we infer that all /j are measurable functions on 
M. By [1], Chapter 2, Theorem 8, any additive measurable function on the 
real line is a continuous linear functional. Hence for each i, we infer the 
existence of Cj G M such that fi{xi) = CiXi holds. Since f{x) = it 
follows that /(x) = (c, x) for some c G Sd- □ 

Also, we consider Cauchy's exponential equation for h : Sj" — )• M+, that is, 

(6.5) h{x + y) = h{x)h{y), x,yeSj. 



Lemma 6.2. Suppose /iiSj"— t-M^. is measurable, strictly positive, and 
satisfies (6.5). Then h{x) = e~'^'^'^\ for some c^Sd- Ifh<l, then cG S^ . 

Proof. Since h is strictly positive, its logarithm yields the well defined 
function g.S^ — >■ M, g{x) :=log/i(x). Clearly g is additive, hence by the 
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first part of Lemma 6.1, there exists a unique additive extension / : S'd — )■ R. 
Also, / is measurable on Sj", hence by the second assertion of Lemma 6.1 
we have f{x) = —{c,x), for some c € Sd- The last statement follows from the 
monotonicity of the exponential and the self duality of Sj". □ 

Remark 6.3. The assumption of strict positivity of h in the preceding 
lemma is essential. Otherwise, there exist solutions h which are not of the 
asserted form. 

Lemma 6.2 is the main ingredient of the proof of the following character- 
ization concerning fc-fold convolutions of Markov processes. 

Lemma 6.4. Let {F^^^)^^g+ eV (i = 0,1, . . . ,k). Then 

(6.6) F%*...*F%=F^^^ Vx»GS+ x = xW + ... + xW, 

if and only if for all t = {h, ...,tN)eR^ andu= {u^^\. . .,u^^^) e (Sj)^ , 
iV G No, there exists < (t, u) < 1 and ij;{t,u) e Sj such that H^Li P^''' i^, 
u) =/>(o)(t,u) and 

(6.7) 

yx£Sj,j = o,i,...,k. 

Proof. We proceed similarly as in the proof of [16], Lemma 10.3. Fix 
A; > 1, > 1, t, u and set 

^(^■)(x):=E(?-)[e-^"i<"^'''^'»>]. 

By the definition of the convolution, (6.6) is equivalent to the following: 

g«(x«) 5('=)(rEW)=gW(x) 

(6.8) 

Hence, the implication (6.7) (6.8) is obvious. For the converse direction, 
we observe that g^'^^ are strictly positive on all of . Thus, by (6.8) we have 

9 := ^(^V^^^) (0) = • • • = 9^'^/9'^'^ (0) = 9^''>/9^'^ (0) 

and (7 is a measurable, strictly positive function on S'^ satisfying (6.5). 
Hence, an application of Lemma 6.2 yields the validity of (6.7), where 
p(»)(t,u) = ^(^^(O). By the definition of g^'\ it follows that < p^'\t,u) < 1 
and V(t,u) □ 

We are prepared to prove Theorem 2.9: 
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(i) =^ (ii): due to Lemma 6.4, infinite decomposability implies that X is 
affine. Also, by the definition of infinite decomposability and by Lemma 
6.4 we have that the fcth root (P^f^) for each k> 1 is an affine process 
with state space with exponents ip{t,u) and (/){t,u)/k. This implies that 

(IPaf^)^g5+ has admissible parameters {a,b/k, /3^^ ,c/k,j,m/k, Hence, the 
admissibility condition proved in Proposition 4.18 implies b/k^ {d— l)a ^ 0, 
for each k, which is impossible, unless a = or d = 1. 

(ii) (iii): follows from Proposition 5.11, in view of the Levy-Khintchine 
form of —4>{t, •) — {ip{t, •),x), for each t> 0. 

(iii) =^ (i): by definition, every transition kernel pt{x, d£,) of X is infinitely 
divisible with Laplace transform P^e"^'*'^^ = e~'^*^*'")~^^''^*^*'")^ For each k > 
1, the maps (p^^^ := ip^^'^ := satisfy the properties (3.1) and (3.2). Also, 
infinite divisibility implies that for each {t,x) G M+ x Sj", 

g(fc)g-(«,x> ._ ^-cl,W{t,u)-{i,(''){t,u),x/k) 

is the Laplace transform of a sub-stochastic measure on 5^. In conjunction 

with Properties (3.1) and (3.2) we may conclude that Q^'^^ gives rise to a 
Feller semigroup on Co(S'j'), which is affine in y = x/k. Hence, we have 
constructed for each fc > 1 a A;th root of X which is stochastic continuous by 
the definition of its characteristic exponents (f)^^\ilj^^\ Thus Theorem 2.9 is 
proved. 

APPENDIX A: EXISTENCE AND VIABILITY OF A CLASS OF 

JUMP-DIFFUSIONS 

In this section, we study existence and viability in a nonempty closed 
convex set D C M" of solutions to the equation 

(A.l) Xt = x+ [ b{Xs)+ [ a{Xs)dWs + Jt, 

Jo Jo 

where b{x) £ Cb(IR"',M"), a{x) £ Cfe(M",]R"^™) are Lipschitz continuous maps, 
W a standard m-dimensional Brownian motion and J a finite activity jump 
process with state-dependent, absolutely continuous compensator K{Xt, d^) dt 
We further assume that x i— )• K{x,M.^) is bounded. 

We tackle this problem in three steps. First, we derive some regularity and 
existence results for diffusion SDEs. These results are not in the standard 
literature, we thus provide full proofs. Second, we prove existence of a cadlag 
solution X for (A.l). Finally, we provide sufficient conditions for X to be 
D- valued. 
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A.l. DifTusion stochastic differential equations. Let J-", (J-f),P) be a 
filtered probability space satisfying the usual conditions and carrying an 
m-dimensional standard Brownian motion W. We consider the following 
diffusion SDE: 

(A.2) Xt = x+ /*6(X,)l|e<4ds+ f a{Xs)l{e<s}dWs, 

Jo Jo 

where {6,x) G [0,oo] x M" and b and a are as above. Recall that X is a 
solution of (A.2) if X is continuous and (A.2) holds for all t > a.s. In 
particular, note that this null set depends on {0,x). 

Lemma A.l. Fix T > and let p>2. Furthermore, let ©1,02 be stop- 
ping times and for i = 1,2, Ui, -measurable random variables. Consider 
the following equations: 

Xt = Ui+ [ b{Xs)l{Q,<s}ds+ [ a{Xs)l{e^<s}dWs, 
Jo Jo 

Yt = U2+ [\{Ys)l{e,<s}ds+ fa{Ys)l{e,<s}dWs. 
Jo Jo 

Then there exists a constant C depending only on p, T, n, the Lipschitz 

constants of h and a and ||6||ooi ||<7||oo such that for 0<t<T, 



E 



snp\\Xs-Ysr 

■s<t 



(A.S) 



< CE 



\\Ui-U2\\p + \@i/\t-e2/\t\p/^ 

+ / supWXu-YuW^ds 

Jo u<s 



Proof. By the same arguments as in the proof of [45], Lemma 11.5, we 
first obtain the following estimate: 

snp\\Xs-Ys\\P 

s<t 

< 3^"' (\\Ui - + \\b{Xs)l{e,<s} - b{Ys)l{e,<s}\\ dsj 

[ {a{Xu)l{e,<u} - f^(^«)l{02<«}) dWu 
Jo 



+ sup 

s<t 



Moreover, 

ft 







\\b{Xs)l{ei<s} - b{Ys)l{e2<s} II ds 
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< 2^' 



(eive2)At \P 

\ms)\\ds) + 

GiAt / 



+ 



(eive2)At \ P 

\\b{Ys)\\ds 

e2At 



\\b{X,)-b(Ys)\\d.s 



(0ive2)At 

\\b{Xs)-b{Ys)rds 



< 2^-1 (^K\ei A i - 62 A t|P + tP^^ 
<K(tP/^\eiAt-e2At\P/^+ [ sup\\Xu-Yu\fds]. 

\ Jo u<s J 

For the stochastic integral part, we apply the Burkholder-Davis-Gundy in- 
equality 

3 PI 

{a{Xu)l{ei<u} - o-(^«)l{02<«}) dWu 



E 



sup 

.s<t 



< KE 



< KE 



lk(^«)l{ei<«} - <y{yu)'^{e2<u}\? duj 



{eive2)At \p/2 
\\a{Xu)fdu] + 

eiAi / 



{eive2)At \p/2 
\\a{Yu)fdu] 

e2Ai 



+ 



\ p/2-| 

\a{Xu) - a{Y^)f du] 



(0ive2)At 



< KE 



<KE 



IGi Ai-e2 At|f/2+ [ \\a{Xs) - a{YsW ds 
Jo 

IGi Af-Os At|f/2+ 1 sn^\\Xu-Yufds 

Jo u<s 



where K always denotes a constant which varies from line to line. The last 
estimate in both inequalities follows from the the Lipschitz continuity of b 
and a. By assembling these pieces, the proof is complete. □ 

Here is a fundamental existence result, which is not stated in this general 
form in the standard literature. Therefore, we provide a full proof. 

Theorem A. 2. There exists a function Z : [0, oo] x M" x O x R+ M" 
with the following properties: 

(i) Z{6,x,uj,t) is continuous in {9,x,t) for alloj. 

(ii) Z is i3([0,oo] x W^) ®V -measurable P 



Here, V denotes the predictable cr-field. 
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(iii) Z{0,x,uj,t) solves ( A. 2) for all {9, x). 

(iv) Let Q be a stopping time and U an Tq, measurable random variable, 
then Xt = Z(Q, U, t) solves 

(A.4) Xt = U+ [\{Xs)l{e<s}ds+ f a{Xs)l{Q<s}dWs. 
Jo Jo 

Proof. For every {0,x) G [0,oo] x M", there exists a unique solution 
Xt(oj) = Z{6,x,uj,t) of (A. 2), which is continuous in t. This is a conse- 
quence of the Lipschitz continuity of x i— )■ 6(x)l|e<s} and x i— )• a{x)l^Q<^gj. 
Uniqueness is meant modulo indistinguishability. From estimate (A. 3), we 
can deduce for p>2,x,ye [-T, T]", < 61,62 <T and <t <T, 



E 



sup\\Z{6i,x, s) 

-s<t 



Z{62,y,s) 



< K 



+ 



E 



sup||Z(ei,x,n)-Z(02,2/,n)|r 



10 ^u<s 

for some constant K. Hence, by Gronwall's lemma, 



ds 



E 



sup||Z(0i,x,s)-Z(e2,2/,s)ir 

s<t 



<Ke^^(||x-yr/2 + |0i-^2r/') 

<cm,x)-{e2,yW'\ 



Let now Dya = {j2 G Z, A; € N} be the set of dyadic rational numbers 
in R and Dya" = Dya x • • • x Dya the set of dyadic rational numbers in M". 
Furthermore, we define M by M = Dya"+^ n([0, T] x [-T,Tf). By setting 
p = 2n + 4, we can apply Kolmogorov's lemma. Indeed, analogous to the 
proof of [32], Theorem 2.8, we derive for all {6i,x), {62,y) G M with < 
11(^1, x) — (02,y)|| < h{u}), where /i is a positive valued random variable, and 
for all uj G fi^, where G is some set depending on T with P(O^) = 1, 
the following estimate: 



(A.5) snp\\Z{0i,x,uj,s) 

s<t 



Zie2,y,u;,s)\\<6\\{di,x)-ie2,yW 



Here, 7 G (0, -) and 5 is some positive constant. Let us now define Z: 
ii oj ^ Q^, then Z{9,x,u},t) = x for < t < T. For u G 0^ and {6,x) G 
M, Z{e,x,u,t) = Z{e,x,u,t) for 0<t<T. If {9,x) G M", we choose a 
sequence {9n,Xn)neN ^ such that {6n,Xn) — )• {d,x). By estimate (A.5), 
Z{dn,Xn,io,t) is a Cauchy sequence converging with respect to supg<;j|| • ||. 
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We can therefore set Z{9,x,uj,t) = lim.n^oo Z{9n,Xn,u},t). As we have uni- 
form convergence in t and as Z is continuous in t, the resulting process Z 
is jointly continuous in {6,x,t). Furthermore, for every {0,x), Z is indistin- 
guishable from Z, that is, 

(A. 6) ¥[Z{e, X, t) = Z{e, X, t) for all < t < T] = 1. 

Indeed, for (6, x) G M, this is clear and for {6, x) G M^, we have for (0„, x„)„,gN ^ 
M with {9n,Xn) {0,x) 

- Zi9,x,s)\\ >e\< Ce-P\\{en,xn) - {0,x)r/^ 

^s<t -I 

which implies that Z{9n.,Xn,t) — )• Z{0,x,t) in probability, uniformly in t. 
As Z{9n.,Xn,t) —7- Z{9,x,t) a.s., and thus in particular in probability, it fol- 
lows that Z{9,x,t) = Z{9,x,t) a.s. for all < f < T. Letting T — t- oo proves 
assertion (i). 

Statement (ii) is then a consequence of (i) and the J-^-measurability of 
u I— 7- Z{9,x,t,oj), which is satisfied since J-q contains all null sets of J-. 

Furthermore, property (A. 6) implies that Z{9,x,t) is a solution of (A. 2) 
for all {9,x), which yields assertion (iii). 

In order to prove (iv), we proceed in two steps: 

Step 1 . We first assume that and U take finitely many values 9i, . . . ,9k G 
[0, oo] and xi, . . . , G M", respectively. Denote 

Aj = {e = 9j}, Bh = {U = Xh}. 

Then 

ZiQ,U,t) = Y,U,nB,Zi9,,Xh,t) 

does the job. Indeed, as Aj n Bh are disjoint and Aj Ci Bh € Tq^ for all j, /i, 
we have (see, e.g., [36], page 39) 

C/+ [\{Z{e,U,s))l{e<s}ds+ f a{Z{e,U,s))l{e<s}dWs 
Jo Jo 

+ f '^'^A,nBhCr{Z{9j,Xh,s))l{g^<sydWs 

j,h 

= '^'^AjnBh(xh+ I b{Z{9j,Xh,s))l{g^<syds 
j,h ^ 
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= J2'^A,nB^z{ej,xh,t) = zie,u,t) 

for alH > a.s. 

Step 2. For general B, U, approximate G*-^'^ IQ hy the simple stopping 
times 

Qik) ^ f i2-^ (j - 1)2"'= < Q < j2"^i = l, . . . , A;2^ 
\oo, A;<Q. 

Let be a sequence of /"e-measurable random variables, each C/^'^ tak- 
ing finitely many values, and U^^^ — )■ [/ in (such U^^^ obviously exists). 
Moreover, {Q^'^) = 9j} n = Xh} G J'oj for aU j,h (see [32], Chapter 1, 
Problem 2.24). 

By Step 1, each Z{Q^^\U^^^) satisfies the respective SDE. Moreover, from 
estimate (A.S) and Grownwall's lemma we deduce that for any T > 0, there 
exists a constant C such that 

E [supll Z(Q W , , t) - z{e^^'\ , t) f " 

< Ce^^E[||C/(*^) - [/(^-'^ f + |G('=) at- G^'^') a T|]. 

Hence, Z(G('=),;7W) is a Cauchy sequence and thus converging with respect 
to E[supj<2-|| • Ip], for all T > 0, to some continuous process X satisfying 
(A. 4). On the other hand, by the continuity of {9,x) i— )• Z{6,x,t), we know 
that 

Z{e^''\U^^\t)^Z{Q,U,t) 

for all uj and t>0. Again, by continuity of 1 1— Z{@, U, t), we conclude that 
Z{@, U) = X up to indistinguishability, which proves the claim. □ 

A. 2. Existence of jump-diffusions. We now provide a constructive proof 
for the existence of a solution of (A.l) on a specific stochastic basis which 
is defined as follows: 

• (f2, -F, iJ^t)t>o) is a filtered space, where := ili x ^l2, '■= Qt ® T~it and 
F = Q®'H are precisely defined below. Note that we do not have a measure 
on {^iJ-) for the moment. The generic sample element will be denoted by 
oj = (wi, W2) S ^■ 

• (r^i, ^, (^t)j>o, Pi) is some filtered probability space satisfying the usual 
conditions and carrying an m-dimensional standard Brownian motion W . 
We shall consider the above diffusion SDE (A. 2) on Vti and thus obtain 
the respective solution Z{9,x,uji,t) satisfying the regularity properties of 
Theorem A. 2. 
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• {il.2,'H) is the canonical space for M"-valued marked point processes (see, 
e.g., [30]): consists of all cadlag, piecewise constant functions a;2 : [0, 
7'oo(w2)) with a;2(0) = and Too{uj2) = lining 00 7'n(w2) < 00, where 

Tn{u}2), defined by Tq = and 

Tni^2) :=inf{t > T„_i(w2) | W2(i) / W2(t-)} A 00, n > 1, 

are the successive jump times of 102- We denote by 



the canonical jump process, and let Tit = (^{Js \ s <t) be its natural fil- 
tration with % = 'Hoo- Note that T„ are (Tit) and (J-()-stopping times if 
interpreted as r„(a;) = r„(a;2). 

The following statement is meant to be pointwise, referring to the filtered 
measure space (Tt)) without reference to a probability measure. 

Lemma A. 3. Let Z{9,x,u}i,t) be as of Theorem A. 2. Then for an J-Tn- 
measurable random variable U{uji,uj2) the process Z{Tn{uj2),U{uJi,u}2),uJi,t) 
is: 

(i) continuous in t for all {uji,u}2), 

(ii) Ft-adapted on {Tn < t}. 

Proof. The first assertion is a consequence of Theorem A.2(i). The 
second one follows from the B{[0, 00] x M"'){g)P-measur ability of Z{9, x,u}i,t), 
as stated in Theorem A.2(ii), and the fact that T„ and U are J^j-measurable 
on {Tn < t}. □ 

Here is our existence result for (A.l). 

Theorem A. 4. There exists a cddldg Tt-adapted process X and a prob- 
ability measure P on {i^,J-) with ¥\g = Pi, such that X is a solution of (A.l) 



Proof. We follow the arguments in the proof of [21], Theorem 5.1, 
which is based on [30], Theorem 3.6, and proceed in three steps. 

Step 1. We start by solving (A.l) along every path uj2- To this end, let us 
define recursively: Aa;2(0) = Aw2(oo) = 0, X^^^ = x, and for n > 1: 



Jt{uj) = Jt{0O2) =^2{t) 



on [0,roo(a;2)) 



on {n,T, {Tt),F). 




t = 00 
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where xq is any fixed point in D C M" and Z satisfies the properties of Theo- 
rem A. 2. By Lemma A. 3, every is continuous in t for all {uji,uj2) and Tt- 

adapted on {Tn < t} since x!p^_^\uji,uj2) + Aa;2(r„_i) is -measurable. 
Thus, the process 



n>l 



is cadlag .Fj-adapted and solves (A.l) on {fli,Q, {Qt),Fi) for t G [0,roo(u;2)) 
and any fixed path uj2- 

Step 2. It remains to show that there exists a probability measure P such 
that K{Xt^d^) is the compensator of J and P|g = Pi holds true. For this 
purpose, we shall make use of [30], Theorem 3.6. Let us define the following 
random measure u by 



v{dt, di) 



(K{Xt,dOdt, t<T, 
1 0, t>T, 



OO ) 



Observe that v is predictable, since Xt is cadlag and J-^-adapted. Theorem 
3.6 in [30] now implies that there exists a unique probability kernel P2 from 
0,1 toTi, such that v is the compensator of the random measure fi associated 
to the jumps of J. On (fi, we then define the probability measure P by 
F{dtj) = Pi((ia;i)P2(ct;i, (ia;2) whose restriction to Q is equal to Pi. 

Step 3. We finally show that X defined by (A.7) solves (A.l) on (0,7", {J^t),l 
for ah t > 0. Note that W{u;) = W{ui) is an (17, J", (Ji),P)-Brownian mo- 
tion. This implies that Z(9,x,uj,t) = Z(6,x,uJi,t) is a solution of (A. 2) on 
(17, (J-t),P), satisfying the properties of Theorem A. 2. It thus remains to 
show that Too = 00 P-a.s. Let fi be the random measure associated with the 
jumps. As X I—)- K{x,W^) is bounded, we have for all T > 0, 



Ep[//([0,r] xM")] =Ep[i/([0,r] xM")]=Ep' 



[ K{Xt,W'^)dt 
Jo 



<CT 



for some constant C. This implies that ^([0,r] x M") < 00 a.s. for all T > 
and hence P[roo < 00] = or equivalently Too = 00 a-s. □ 



A.S. Viability of jump-diflFusions. Consider a nonempty closed convex 
set D C M". We now provide sufficient conditions for the solution X in 
(A.7) to be D-valued. This result is based on [13], Theorem 4.1. We recall 
the notion of the normal cone 

(A.S) Nd{x) = {u£ M" \{u,y-x)> 0, for all y e D} 

of D at X ^ D, consisting of inward pointing vectors. See, for example, [28], 
Definition III. 5. 2. 3, except for a change of the sign. 
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Theorem A. 5. Assume that a also has a Lipschitz continuous deriva- 
tive. Suppose furthermore that 

(A.9) x + supp{K{x,-))CD, 

(A.IO) {a'{x),u)=0, 

(A.ll) (^b{x)-^f2Da\x)a\x),u^>0, 

for all u G No{x) and x & D, where denotes the ith column of a. Then, 
for every initial point x G D, the process X defined in (A. 7) is a D-valued 
solution of (A.l). 



Proof. We have to show that Xt = X]„>i X^ '^(a;i,c<J2)l{T„_i<t<T„} £ D 

a.s. for all t > 0. We proceed by induction on n. For n = 1, x[^\ is simply 
given by 

Xf)=x+ fb{x'^P)ds+ f a{x'^P)dWs. 
Jo Jo 

Due to [13], Theorem 4.1, conditions (A.IO) and (A.ll) imply that for all 

t >0, XI £ D a.s. Let us now assume that for all t > 0, £ D a.s., 

thus in particular X^ = Xt„_^- E D a.s. If T^-i = oo, then we immedi- 
ately obtain 

Otherwise, let / G Cb(M",M+) satisfy supp(/) C D". Then, 



E 



/ f{XT,,_,^+OK[XT,,_,^,dO 

Jk."\{o\ 



0, 



/R"\{0} 

since by (A.9), Xt^_^- + suy>v{K {Xt„_^- , ■)) ^ D a.s. and f{D) = 0. Hence, 
/(^ti^ + ^JTn-i) = a.s., implying that X^^'^^ + AJt„_i ^ supp(/) a.s. 
As this holds true for all / G C6(M",]R+) with supp(/) C D^, it follows that 
X^^Zl'^ + A Jt„_i G D a.s. Thus, again by [13], Theorem 4.1, and conditions 
(A.IO) and (A.ll) 

4)=47f + AJr„_,+ A(4"))l|^_,<,|d.+ fa{X(-%T„^^<,ydWs 

Jo Jo 

a.s. takes values in D, which proves the induction hypothesis. The definition 

of X then yields the assertion. □ 
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APPENDIX B: AN APPROXIMATION LEMMA ON THE CONE OF 

POSITIVE SEMIDEFINITE MATRICES 

In this section, we deliver a difFerentiable variant of the Stone- Weierstrass 
theorem for C°°-functions on Sj". This approximation statement is essential 
for the description of the generator of an affine semigroup, as is elaborated 
in Section 4.2. 

We employ multi-index notation in the sequel. For n > 1, a multi-index 
is an element o; = (ai, . . . , a„) G Nq having length \a\ := ai -|- • • • -|- a„. The 
factorial is defined by a! := HiLi'^*" "^^^ partial order < is understood 
componentwise, and so are the elementary operations -|-, — . That is, a> (3 
if and only if Oi > Pi for i = 1, . . . , n. Moreover, for a> /3, the multinomial 
coefficient is defined by 

f a\ _ al 
[f3) - (a-/3)!/3!- 

We define the monomial := Yl^^i and the differential operator 9" := 
af la„ ■ Corresponding to a polynomial P{x) = ^|q,|</c a^x", we introduce 

the differential operator P{d) := J2\a\<k'^ct^^- 

Let S = S{Sd) denote the locally convex space of rapidly decreasing C°°- 
functions on (see [46], Chapter 7), and define the space of rapidly de- 
creasing C°° -functions on 5^" via the restriction 

(B.l) S+ = {f = F\,y.FGS}. 

Equipped with the increasing family of semi-norms 

(B.2) Pk,+ {f)-= sup |x"a^/(x)|, 

xeS+ ,\a+(}\<k 

5+ becomes a locally convex vector space (see [46], Theorem 1.37). 
For technical reasons, we also introduce for e > the semi-norms 

PkM)-= sup |x"a^/(x)| 

xf^S++B<,{0),\a+l3\<k 

on C°°{S(i), where B<r{y) = {z £ Sd \ \\z — y\\ < r} denotes the closed ball 
with radius r and center y. Note that pk,+ =Pk,o- We first give an alternative 
description of 5+. 

Lemma B.l. We have 
S+ = {f = G\g+:Ge C^iSd) and 3e > such that Pk,e{G) < oo VA; > 0}. 
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Proof. The inclusion C is trivial. Hence, we prove 5. So let f = G \ q+ 

for some G £ G°°{Sd) with pk,e{G) < oo for all A; > and some e > 0. 

We choose a standard mollifier p G C^{Sd) supported in -B<e/2(0) and 
satisfying p > 0, / /> = 1- For 6 > we introduce the neighborhoods Kg := 
S'^ + B<s{0) of Sj". The convolution 93 := Ij^^^^ GC°°(S'd) of the indicator 

function for with p satisfies 93 = 1 on Sj" and it vanishes outside K^. 
Furthermore, all derivatives of (p are bounded, since 



d"p{y - x) dy 



d"p{z) dz 



< / \d'^p{z)\dz < 00, 

where the last estimate holds because supp/O C i?<e/2(0)- 

Now we set F :=G -(p. By construction F G C°°{Sd), F\g+= f and F 

vanishes outside K^, because ip does. What is left to show is that F £ S. 

Since F vanishes outside K^, it is sufficient to deliver all estimates of its 

derivatives on Kg;. 

Let a,P £ then we have by the Leibniz rule 



= E (^){d''-''y^{x)){x'^d''G{x)). 



0<7</3 ^ ' 

By assumption d'^G is bounded on K^, and {d^~'^ ip{x)) is bounded on all 
of Sd- Hence, by the last equation, we have sup^.g^^ |q,^^|<;^|x" dl^F{x)\ <oo 
for all /c G No, which by definition means F £ S. □ 

Lemma B.2. Let u G S'^^ . Then for each £ >0, and for all k>0 we 
have pfc^e(exp(— (t/, •))) <oo. In particular, we have 

fu :=exp(-(n, g5+. 
That is, fu = Fu\ <?+ for some Fu £S. 

Proof. Since u G -Sj"*", there exists a positive constant c such that 
{u,x) > c\\x\\, for all x G S'j'. Hence, we obtain by a straightforward cal- 
culation, pfc^+(exp(— (n, •))) < 00, for all A: > 0. 

Next, let e > 0, and write x = y + z, where y G Sj" and z G i?<e(0) and 

pick multi-indices q;,/3 G J^jj^'^"'''^^/^. Then we have by the binomial formula 



X 
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0<7<q; ^'^^ 

Now since z ranges in a compact set, and since pfc^4.(exp(— (n, •))) < oo we 
see that x°^d^e~^'^'^'^ must be bounded uniformly in x G +i?<e(0). Hence, 
Pfc^£(exp(— (ti, •))) < oo, for all A; > 0. 

Together with Lemma B.l, this implies fu G 5+. □ 

We are now prepared to deliver the following density result for the M- 
linear hull of {/„ = exp(— (n, •))|^+ , u G Sj""*"} in 5+. 

Theorem B.3. M is dense in S+. 

Proof. Denote by 5' = S'{Sd) and 5^ the topological dual of S and 
5+, respectively. The former, 5', is known as the space of tempered distri- 
butions. The distributional action is denoted by (•,•) and (•,•)+ for S' and 
S'^, respectively. 

Now suppose by contradiction, that 7W is not dense in 5+. Then by [46], 
Theorem 3.5, there exists some Ti G S'_^ \ {0} such that Ti = on Hence, 
{Ti, fu)+ = 0, for all u G The restriction F i— F\g+ yields a continuous 
linear embedding S 5+. Hence, the restriction T of Ti to 5, given by 

{T,(p) := {Ti,Lp\g+) + , LpeS{Sd), 

yields an element of S' with supp(r) C Sj". Pick an G 5 according to Lem- 
ma B.2. By the definition of T, we have (T, F„) = {Ti,fu)+ = 0, for all u G 
Sj"*". By the Bros-Epstein-Glaser theorem (see [42], Theorem IX. 15), there 
exists a function G G G{Sd) with supp(G) C Sj", polynomially bounded [i.e., 
for suitable constants C,N we have |G(a;)| < C(l + ||x||)^, for all x G 5j"] 
and a real polynomial P{x) such that P[d)G = T in S' . Hence, we obtain 
for any u G S'j""'' 

= (T, Fu) = {P{d)G, Fu) = {G, P{-d)Fu) 

= [ G{x)P{-d)Fu{x)dx = P{u) [ G{x)exp{-{u,x))dx. 
Js+ Js+ 

But the last factor is just the Laplace transform of G. This implies G = 0, 
hence T = 0, which in turn implies that Ti vanishes on all of 5+, a contra- 
diction. □ 

Acknowledgments. We thank Martin Keller-Ressel and Alexander Smirnov 
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